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1. INTRODUCTION 

The editors suggested a review entitled "Are There Large-Scale Motions in the 
Universe?" . The answer is "yes" , in the sense that the interpretation of the data as 
motions is the simplest model, so far consistent with all other available data under the 
current "standard model" of physical cosmology. I will review tests which could have ended 
up falsifying this model and failed, but the scope of this review is much extended as the 
field has developed far beyond the question of existence of motions. Having adopted the 
motions as a working hypothesis, the study of large-scale dynamics is becoming a mature 
scientific field where observation and theory are confronted in a quantitative way. It is this 
area of major activity in cosmology that is addressed here. 

I have made no attempt to provide a complete reference list, nor have I tried a balanced 
discussion of all the issues of relevance and all the authors involved. My goal is to provide a 
critical account of some of the issues which I find important, with emphasis on theoretical 
implications. In many cases I am not careful in giving proper credit by quoting only a 
recent paper automatically implying "and references therein" . The reader is referred to a 
comprehensive, observation-oriented review of large-scale motions in historical perspective 
by Burstein (1990b), a detailed review of distance indicators in a collection of essays by 
Jacoby et al. (1992), and to "Principles of Physical Cosmology" by Peebles (1993). I 
apologize for tending to describe in more detail work which I was involved in and therefore 
know better. 

The current phase of the field was seeded by two major developments. One was the 
confirmation of the dipole moment in the Cosmic Microwave Background (CMB) (Corey 
& Wilkinson 1976; Smoot et al. 1977), indicating via Doppler shift that the Local Group 
of galaxies (LG) is moving at ~ 600 kms~"^ relative to the cosmological frame defined by 
the CMB. The other was the invention of methods for inferring distances independent of 
redshifts based on intrinsic relations between galaxy quantities (§3; TuUy & Fisher 1977, 
TF; Faber & Jackson 1976, FJ). The radial peculiar velocity of a galaxy (the "velocity" u) 
is the difference between its total radial velocity as read from the redshift (the "redshift" z) 
and the Hubble velocity at its true distance (the "distance" r). Improved versions of these 
methods reduced the distance errors to the level of 15-21% which, with several hundred 
measured galaxies across the sky, enabled modeling the large-scale velocity field in terms 
of few-parameter "toy" models (§4.1), starting with a Virgo-centric infall (Aaronson et 
al. 1982b) and ending with spherical infall into a "Great Attractor" (GA, Lynden-Bell 
et al. 1988). The finding by the "seven samurai" (7S, Burstein et al. 1986) that the LG 
participates in a large streaming motion launched the present high-profile activity in this 
field. The toy modeling is gradually being replaced by non-parametric methods, where the 
full velocity field is reconstructed based on properties of gravitational flows (§4) and the 
associated mass-density fluctuation fleld is recovered from the spatial velocity derivatives 
(§2). With no simplifled geometry imposed, the motions are not associated with single 
speciflc "sources" ; the gravitational acceleration is an integral of a continuous density fleld 
consisting of swells and troughs simultaneously pulling and pushing. 

A parallel major development has been of all-sky magnitude-limited redshift surveys 
with many thousands of galaxies, starting with the CfA and SSRS optical surveys and 
continuing with the very useful recent surveys based on the IRAS satellite (§5). The 



large-scale inhomogeneity in the galaxy distribution {e.g. de Lapparent et al. 1986) 
provided a clear hint for associated motions. An all-sky redshift survey can be converted 
into a galaxy-density field and then integrated to derive a predicted velocity field under 
the assumption of gravity and a certain "biasing" relation between galaxies and mass. 
The comparison of the fields obtained from redshifts to those obtained from velocities is 
at the heart of the research of large-scale structure (LSS), and the results carry major 
implications (§6.2, 8.2). 

Data of both types are rapidly accumulating, and a major effort is directed at reducing 
the errors and carefully estimating those which remain, to enable quantitative testing 
of LSS formation theories. The standard theory consists of several working hypotheses 
which one tries to falsify by the observations or, if found consistent, to determine the 
characteristic model parameters. The hypotheses, which will be elaborated on later, can 
be listed as follows: 

HI. The background cosmology is the standard homogeneous Friedman Robertson Walker 
model, possibly with an Infiation phase, where the CMB defines a cosmological "rest 
frame" . If so, then one wishes to determine the cosmological density parameter O 
(and the cosmological constant A and the Hubble constant H). 

H2. The structure originated from a random field of small-amplitude initial density 
fluctuations. If so, the goal is to find out whether they were Gaussian, whether 
the power spectrum (PS) was scale-invariant (power index n = 1), and whether the 
energy density was perturbed adiabatically or in an isocurvature manner. 

H3. The spectrum of fiuctuations was filtered during the radiation-plasma era in a way 
characteristic of the nature of the dark matter (DM) which dominates the mass density. 
The DM could be baryonic or non-baryonic. If non-baryonic it could be "hot" or "cold" 
depending on when it became non-relativistic. 

H4. The fiuctuations grew by gravitational instability (GI) into the present LSS. This is a 
sufficient but not necessary condition for: 

H4a. The quasi-linear velocity field smoothed over a sufficiently-large scale is 
irrotational. 

H4b. The galaxies trace a unique underlying velocity field, apart from possible "velocity 
bias" of ~ 10% on small scales. 

H5. The density fiuctuations of visible galaxies are correlated with the underlying mass 
fiuctuations. If this relation is roughly linear, then the linearized continuity equation 
in GI implies a relation between velocity and galaxy density. If so, the characteristic 
parameter is the density-biasing factor b. 

H6. The TF and Dn—cr methods measure true distances^ which allow the reconstruction 
of a large-scale velocity field with known systematic biases under control. 

This review is geared toward the confrontation of observations with these hypotheses. 
The relevant observations can be classified into the following three major categories: 

01. Angular fiuctuations in the CMB temperature at various angular scales. 

02. The distribution of luminous objects on the sky and in redshift space. 
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03. Peculiar velocities of galaxies along the line of sight. 

Note that 01 and 03 are related to dynamical ingredients Hl-4 and H6, bypassing the 
uncertain nature of galaxy-density biasing H5. Also, 01 and 02 refer to the theory 
independently of H4a,b and H6, which address the velocities and their analysis. 

2. GRAVITATIONAL INSTABILITY 

This section provides a brief account of the standard theory of GI, and of the linear and 
quasi-linear approximations which serve the analysis of motions. Let s,«, and $(, be the 
position, peculiar velocity and peculiar gravitational potential in comoving distance units, 
corresponding to ax, av, and a^^g in physical units, with a{t) the universal expansion 
factor. Let the mass-density fluctuation he 5 = {p — p)/p- The equations governing 
the evolution of fluctuations of a pressureless gravitating fluid in a standard cosmological 
background during the matter era are the Continuity equation, the Euler equation of 
motion, and the Poisson fleld equation {e.g. Peebles 1993): 



(^ + V-w + V-(w(^) = , (1) 

v + 2Hv + {v-V)v = -V^g , (2) 

V^$g = (3/2)ij2 0(^ , (3) 



with H and O varying in time. The dynamics do not depend on the value of the Hubble 
constant H; it is set to unity by measuring distances in kms~"^ (1 h~"^Mpc = 100 kms~"^). 

In the linear approximation, the GI equations can be combined into a time evolution 
equation, 8 -\- 2H5 = {3/2)H^Q5. The growing mode of the solution, D{t), is irrotational 
and can be expressed in terms of /(O) = H~^ D/D sa 0°'^ (see Peebles 1993, eq. 5.120). 
The linear relation between density and velocity is 

S = So = -{Hf)-'V-v . (4) 

The use of 5q is limited to the small dynamical range between a few tens of megaparsecs 
and the ~ 100 h~"^Mpc extent of the current samples. In contrast the sampling of galaxies 
enables reliable dynamical analysis with smoothing as small as ~ 10 h~"^Mpc, where |V-w| 
obtains values > 1 so quasi-linear effects play a role. Unlike the strong non-linear effects 
in virialized systems which erase any memory of the initial conditions, mild non-linear 
effects carry crucial information about the formation of LSS, and should therefore be 
treated carefully. Figure 1 shows that 5q becomes a severe underestimate at large \5\. This 
explains why eq. (4) is invalid in the non-linear epoch even where 5 = 0; the requirements 
that J5d^x = by deflnition and J^-vd^x = by isotropy imply — V-«>(^ at \5\ ^ 1. 
Fortunately, the small variance of V-w given 5 promises that some function of the velocity 
derivatives may be a good local approximation to 5. 

A basis for useful quasi-linear relations is provided by the Zel'dovich (1970) 
approximation (Z). The displacements of particles from their initial, Lagrangian positions 
q to their Eulerian positions x at time t are assumed to have a universal time dependence, 

x{q,t) - q = D{t)i^{q) = f-'v{q,t) . (5) 
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Figure 1: Quasi-linear velocity-to-density approximations. A = Sapproxi^) ~ 
Strue- The mean and standard deviation are from large standard-CDM N-body 
simulations normalized to erg = 1, Gaussian-smoothed with radius 12 h~"^Mpc 
(see Mancinelli et al. 1994). Note the factor of 5 difference in scale between the 
axes. 

For the purpose of approximating GI, the Lagrangian Z approximation can be interpreted 
in Eulerian space, 

q{x) = X — f~^v{x) , (6) 

provided that the flow is laminar [i.e. that multi-streams are appropriately smoothed over). 
The solution of the continuity equation then yields (Nusser et al. 1991) 



Sc{x) = \\I — f -^dv/dxW — 1 



(7) 



where the bars denote the Jacobian determinant, and / is the unit matrix. The Z 
displacement is flrst order in f~^ and v, so 5^ involves second- and third-order terms 
(?7Zi,2, ?^u3) as well. The relation (7) is not easily invertible to provide V-w or v when 5 is 
given, but a useful approximation derived from simulations is 



V-w 



-fS/{l + 0.18S) . 



A modifled approximation, which is derived by adding a second-order term to the Z 
displacement (Moutarde et al. 1991) and truncating all the expressions at second order 
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while solving the continuity equation is (Gramann 1993a) 



>c2 



-f ^V-w + (4/7)/ ^m„2 , m^2 = ^^{^^v,^JVJ-^JV,^,VJ). (9) 

i j>i 



The factor 4/7 replaces 1 in the second-order term of Sc- While terms are kept to 
second-order, it is still not an exact solution to the second-order equations of GI. This 
relation can be inverted in second-order to provide V-t> given 5^ with 112^2 replaced by an 
analogous expression mg2^ involving the gravitational acceleration g. 

Since the variance of 5 given V • t> is small, one expects that a non-linear function 
of V-« which properly corrects for the systematic deviation can be a good quasi-linear 
approximation to 5 (and vice versa). Assuming Gaussian initial fluctuations, Bernardeau 
(1992) found a solution in the limit of vanishing variance: 

S, = [l-{2/3)f-'V-v)]'/'-l, (10) 

which is easily invertible. A polynomial expansion with non-vanishing variance should 
have the form (Zehavi & Dekel, in prep.) 

S^{V.v) = -f-'V-v + a2f-'[{V-vf-IJ2]+asf-'[{V-vf-iJs] + ... . (11) 

Because the flrst two terms vanish when integrated over a large volume, the moments 
jjn = ((V-w)") must be subtracted off to make the n* -order term vanish as well. The 
coeflicients can be crudely approximated analytically {e.g. Bernardeau 1992) or, using 
CDM simulations and Gaussian smoothing, the best coeflicients are 02 ~0.3 and as sa —0.1, 
tested for O values 0.1 — 1 and smoothing radii 5 — 12 h~"^Mpc at erg = 1 (erg is the rms 
of unsmoothed mass-(^ in top-hat spheres of radius 8 h~"^Mpc). The structure of eq. 
(11) makes it robust to uncertain features such as O, the shape of the fluctuation power 
spectrum, and the degree of non-linearity as determined by the fluctuation amplitude and 
the smoothing. Such robustness is crucial when using a quasilinear approximation for 
determining O, for example (§8). 

Fig. 1 demonstrates the accuracy of the explicit quasi-linear approximations using 
CDM N-body simulations and 12 h~"^Mpc smoothing (Mancinelli et al. 1994). 8c^ of 
scatter ~0.1, is an excellent approximation for 8>1 but it is a slight overestimate at the 
negative tail. 8(.2 and 8}, do better at the negative tail, but they are severe underestimates 
in the positive tail. 8s{^ -v) is an excellent robust flt over the whole quasi-linear regime. 
(^c3 is constructed from the three terms in the expansion of 8c in powers of f~^ but with 
the numerical coeflicients adjusted to achieve best flt in the simulation ( — 1.05,0.9,1.5 
replacing unity, independent of O). 



3. MEASURING PECULIAR VELOCITIES 
3.1. Distance Indicators 

Measuring redshift-independent distances to many galaxies at large distances is the key 
to large-scale dynamics (review: Jacoby et al. 1992). The simplest method assumes that 
a certain class of objects is a "standard candle" , in the sense that a distance-dependent 
observable is distributed intrinsically at random with small variance about a universal 
mean. The luminosity of an object (oc r~^) or its diameter ((oc r~"^), can serve as this 
quantity. In a pioneering study, Rubin et al. ( 1976a, b) used the brightness of giant Sc 
spirals to discover a net motion for the shell at 35—60 h~"^Mpc that agrees within the errors 
with more modern results, but the large uncertainties in this simple distance indicator made 
this result controversial at the time. 

So far, the most useful distance indicators for LSS have been of the TF-kind, based 
on intrinsic relations between two quantities: a distance-dependent quantity such as the 
flux (xLlr'^^ and a distance-independent quantity a - the maximum rotation velocity of 
spirals (TF) or the velocity dispersion in ellipticals (FJ). The intrinsic relations are power 
laws, Loca^ , i.e. 

M{r]) = a-br] , (12) 

where M = — 2.51ogL -|- const is the absolute magnitude and r] = logcr. The slope b can 
be determined empirically in clusters, where all the galaxies are assumed to be at the 
same distance, typically yielding /3sa3 — 4, depending on the luminosity band {e.g. /3jsa3, 
/3jf ~ 4). Then, for any other galaxy with observed r] and apparent magnitude m = 
— 2.51og(L/r^) -\- const, one can determine a relative distance via 5logr = m — M {rj) . There 
exists a fundamental freedom in determining the zero point, a, which flxes the distances at 
absolute values (in kms~"^, not to be confused with H which translates to Mpc). Changing 
a, i.e. multiplying the distances by a factor (1 + e) while the redshifts are flxed, is equivalent 
to adding a monopole Hubble-like component —er to w, and an offset 3e to 5 (eq. 4). It 
has been arbitrarily determined in several data sets, e.g. by assuming u = for the Coma 
cluster, but a is better determined by minimizing the variance of the recovered peculiar 
velocity fleld in a large "fair" volume. The original TF technique has been improved 
by moving from blue to near-infrared photometry (H band, Aaronson et al. 1979) and 
recently to CCD R and I bands, where spiral galaxies are more transparent and therefore 
the intrinsic scatter is reduced to a^ ~0.33 mag, corresponding to a relative distance error 
of A= (lnl0/5>7^?^0.15. 

A distance indicator of similar quality for ellipticals has proved harder to achieve. 
Minimum variance, corresponding to A = 0.21, was found for a revised FJ relation involving 
three physical quantities: DI" oc a^ with D the diameter and / oc L/D^ the surface 
brightness (Dressier et al. 1987; Djorgovski & Davis 1987). The parameters were found to 
be a sa 5/6 and (3 sa 4/3. By deflning from the photometry a "diameter" at a flxed value 
of enclosed /, termed Dn, the relation returns to a simple form, Dn oc a^ , similar to FJ 
but with reduced variance. 

The physical origin of the scaling relations is not fully understood, reflecting our 
limited understanding of galaxy formation. What matters for the purpose of distance 
measurements is the mean empirical relation and its variance. However, one can point at 



an important physical difference between the two relations (Gunn 1989), which is relevant 
to the testing for environmental effects (§6.3). The D^—fJ relation is naturally explained 
by virial equilibrium, a^ (xM/D, and a smoothly varying M / L (x M~'^ , which together yield 
£)/1/(i+t) oc a^^^ u/C^+i)^ but the TF relation, involving only two of the three quantities 
entering the virial theorem, is more demanding - it requires an additional constraint which 
is probably imposed at galaxy formation. 

There is some hope for reducing the error in the TF method to the ~ 10% 
ball park by certain modifications, e.g. by restricting attention to galaxies of normal 
morphology (Raychaudhury 1994). But the most accurate to date is the estimator based 
on surface-brightness fluctuations (SBF) in ellipticals (Tonry 1991), where the standard 
candle is the luminosity function of bright stars in the old population. These stars show 
up as distance-dependent fluctuations in sensitive surface-brightness measurements. The 
technique is being applied successfully out to ~ 30 h~"^Mpc [e.g. Dressier 1994), with the 
improved accuracy of ~ 5% enabling high-resolution non-linear analysis, and it can be of 
great value for LSS if applied at larger distances. The need to remove sources of unwanted 
fluctuations such as globular clusters requires high resolution observations which could be 
achieved by HST or adaptive optics. 

The prospects for the future can be evaluated by estimating the length scale over 
which LSS dynamics can be studied using a distance indicator of relative error A. The 
error in a velocity derived from N galaxies at a distance ~ r is ay ~ rA/viV. Let the 
mean sampling density be n. Let the desired quantity be the mean velocity V in spheres of 
radius R, and assume that its true rms value is V20 at i? = 20 h~"^Mpc and oc R~^^~^^' on 
larger scales, with n the effective power index of the fluctuation spectrum near R. Then 
the relative error in V is 

where distances are measured in h~"^Mpc. The observations indicate that V20 ~500 kms~"^ 
and n ~ —0.5 for R = 20 — 60 h~"^Mpc (§7.1). Thus, with ideal sampling of n ~ 
0.01 (h~"^Mpc)~^, the relative error is always only a few percent of r I R. This means that 
LSS motions can in principle be meaningfully studied at all distances r with smoothing 
R ~ O.lr, as long as n ~ —0.5 at the desired R. Since n seems to be negative out 
to ~ 100 h~"^Mpc (§7.1), dense deep TF samples are potentially useful out to several 
hundred megaparsecs. However, several technical difficulties pose a serious challenge at 
such distances. For example, the calibration requires faint cluster galaxies which are harder 
to identify, aperture effects become severe, the spectroscopy capability is limited. 

3.2. Malmquist Biases 

The random scatter in the distance estimator is a source of severe systematic biases 
in the inferred distances and peculiar velocities, which are generally termed "Malmquist" 
biases but should carefully be distinguished from each other {e.g. Lynden-Bell et al. 1988; 
Wilhck 1994a,b). 

The calibration of the TF relation is affected by the selection bias (or calibration 
bias). A magnitude limit in the selection of the sample used for calibration at a flxed true 
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distance {e.g. in a cluster) tilts the "forward" TF regression line of M on ry towards bright 
M at small r] values. The bias extends to all values of r] when objects at a large range of 
distances are used for the calibration. This bias is inevitable when the dependent quantity 
is explicitly involved in the selection process, and it occurs to a certain extent even in the 
"inverse" relation rji^M) due to existing dependences of the selection on rj. Fortunately, 
the selection bias can be corrected once the selection function is known {e.g. Willick 1991; 
1994a). 

The TF inferred distance, c/, and the mean peculiar velocity at a given c/, suffer from 
an inferred- distance bias, which we term hereafter "M" bias. I comment later (§4.4) on 
a possible way to avoid the M bias by performing an inverse analysis in 2;-space, at the 
expense of a more complicated procedure and other biases. Here I focus on a statistical 
way for correcting the M bias within the simpler forward TF procedure in (/-space. This 
bias can also be corrected in an inverse TF analysis in (/-space, using the selection function 
S{d) which is in principle derivable from the sample itself (Landy & Szalay 1992). 

The current POTENT procedure uses the forward TF relation in (/-space. If M 
is distributed normally for a given ry, with standard deviation cr^, then the TF-inferred 
distance c/ of a galaxy at a true distance r is distributed log-normally about r, with relative 
error Asa0.46(Tm. Given c/, the expectation value of r is {e.g. Willick 1991): 

C rP{r\d)^r .IT ^''^(O exp (- ^'"^/^f^^' ) dr 

E{r\d) = -^"oo . ' = 7 { , 14 

^ ' j- P{r\d)dr J- ,2^(,) exp (-lii^^^) dr ' 

where n{r) is the number density in the underlying distribution from which galaxies were 
selected (by quantities that do not explicitly depend on r). The deviation of E{r\d) from 
(/ reflects the bias. The homogeneous part (HM) arises from the geometry of space - the 
inferred distance c/ underestimates r because it is more likely to have been scattered by 
errors from r > d than from r < d, the volume being oc r^. If n = const, equation (10) 
reduces to 

E(r|(/) =(/e'-'^' , (15) 

in which the inferred distances are simply multiplied by a factor, 8% for A = 0.15, 
equivalent to changing the zero-point of the TF relation. The HM bias is corrected this 
way on a regular basis since 7S. 

Fluctuations in n{r) are responsible for the inhomogeneous bias (IM), which is worse 
because it systematically enhances the inferred density perturbations and the value of O 
inferred from them. If n{r) is varying slowly with r, and if A ^ 1, then eq. (14) reduces 
to 



E{r\d) =d 



.9 .9 /c/lnn\ 
1 + 3.5A' + A' H^ 
\ cilnr I 



r^d- 



(16) 



showing the dependence on A and the gradients of n{r\ To illustrate, consider a lump of 
galaxies at one point r with u = 0. Their inferred distances are randomly scattered to the 
foreground and background of r. With all galaxies having the same z, the inferred u on 
either side of r mimic a spurious infall towards r, which is interpreted dynamically as a 
spurious over density at r. 
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In the current data for POTENT analysis (§4.2) IM bias is corrected in two steps. 
First, the galaxies are heavily grouped in z-space (Willick et al. 1994), reducing the distance 
error of each group of N members to A/VN and thus significantly weakening the bias. 
Then, the noisy inferred distance of each object, c/, is replaced by E[r\d) (eq. 14), with 
an assumed n[r) properly corrected for grouping. This procedure has been tested using 
realistic mock data from N-body simulations (Kolatt et al. , in prep.), showing that IM 
bias can be reduced to a few percent. The practical uncertainty is in n[r)^ which can be 
approximated by the high-resolution density field of IRAS or optical galaxies (§5), or by 
the recovered mass-density itself in an iterative procedure under some assumption about 
how galaxies trace mass. The second-step correction to 5 recovered by POTENT is <20% 
even at the highest peaks (Dekel et al. 1994). 

3.3. Homogenized Catalogs 

Several samples of galaxies with TF or D^—fJ measurements have accumulated in the 
last decade. Assuming that all galaxies trace the same underlying velocity field (§6.3), 
the analysis of large-scale motions greatly benefits from merging the different samples into 
one self-consistent catalog. The observers differ in their selection procedure, the quantities 
they measure, the method of measurement and the TF calibration techniques, which cause 
systematic errors and make the merger non-trivial. The original merged set, compiled by 
D. Burstein (Mark II) and used in the first application of POTENT (Bertschinger et al. 
1990), consisted of 544 ellipticals and SO's (Lynden-Bell et al. 1988; Faber et al. 1989; 
Lucey & Carter 1988; Dressier & Faber 1991) and 429 spirals (Aaronson et al. 1982a; 
Aaronson et al. 1986; 1989; Bothun et al. 1984). The current merged set (Willick et al. 
1994, Mark III of the Burstein series) consists of ~ 2850 spirals (Mark II plus Han & Mould 
1990, 1992; Mould et al. 1991; Wilhck 1991; Courteau 1992; Mathewson et al. 1992) and 
the ellipticals of Mark II. This sample enables a reasonable recovery of the dynamical fields 
with ~ 12 h~"^Mpc smoothing in a sphere of radius ~60 h~"^Mpc about the LG, extending 
to ~80 h~"^Mpc in certain regions (§4). Part of the data are shown in Figure 2. 

As carried out by this group, merger of catalogs involves the following major steps: (1) 
Standardizing the selection criteria, e.g. rejecting galaxies of high inclination or low r] which 
are suspected of large errors and sharpening any z cutoff. (2) Rederiving a provisional TF 
calibration for each data set using Willick's algorithm (1994a) which simultaneously groups, 
fits and corrects for selection bias, and then verifying that inverse-TF distances to clusters 
are similar to the forward-TF distances. (3) Starting with one data set, adding each new 
set in succession using the galaxies in common to adjust the TF parameters of the new set if 
necessary. (4) Using only one measurement per galaxy even if it was observed by more than 
one observer to ensure well defined errors, and using multiple observations for a "cluster" 
only if the overlap is small {e.g. <50%). (5) Adding the ellipticals from Mark II, allowing 
for a slight zero-point shift (§6.3). Such a careful calibration and merger procedure is 
crucial for reliable results - in several cases it produced TF distances substantially different 
from those quoted by the original authors. 

4. ANALYSIS OF OBSERVED PECULIAR VELOCITIES 

Given radial peculiar velocities Ui sparsely sampled at positions Xi over a large volume, 
with random errors cr^, the first goal is to extract the underlying three-dimensional velocity 
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Figure 2: Inferred radial peculiar velocities of grouped galaxies in a ±20° slice 
about the Supergalactic plane from the homogenized Mark III catalog (Willick et 
al. 1994). Distances and velocities are in 1000 kms~^. The area of each circle 
marking the object position is proportional to the object richness. This slice 
contains 453 objects made of 1124 galaxies out of 1214 objects in the whole volume. 
Solid and dashed lines distinguish between outgoing and incoming objects. The 
positions and velocities are corrected for IM bias. Note the GA convergence (left) 
and the PP convergence (right-bottom). 

field, v{x). Under GI, this velocity field is subject to certain constraints, e.g. being 
associated with a mass-density fluctuation fleld, 8{x) - the other target for recovery. A fleld 
can be deflned by a parametric model or by the fleld values at grid points. The number of 
independent parameters or grid points should be much smaller than the number of data 
points in view of the noisy data. 

4.1. Toy Models 

Given a model velocity fleld t>(ayt,^), the free parameters Q,k can be determined 
globally by minimizing a weighted sum of residuals, e.g. 



-logL oc ^ Wj [u. 



X, ■v{ak,Xi)f 



(17) 



If the errors are Gaussian and Wi oc af then L is the likelihood, and it is a useful 
approximation for log-normal errors as well. The model could flrst be a few-parameter 
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"toy" model with simple geometry. Already the simplest bulk-flow model, v[x) = B 
corresponding to 5 = 0^ is of interest because the data clearly show a bulk flow component 
of several hundred km s~^ in our neighborhood (§7.1). Another simple model is of spherical 
symmetry, expected to be a reasonable flt in voids and in regions dominated by one high 
density peak (Bardeen et al. 1986). The velocity proflle as a function of distance r from 
the infall center is not particularly constrained by GI, and a speciflc proflle which proved 
successful is 



v\n = -vig 






-| (n+l)/2 



The center is specifled by its angular position and its distance rig from the LG, and the 
proflle is characterized by its value vig at the LG, a core radius Tc, and a power index n. 
For r <^ r,. the velocity rises oc r and for r ^ Tc it falls off oc r~". The associated density 
proflle is given in the linear approximation by the divergence 8 = — f~^ r~^ dr[r v{r)\. The 
7S ellipticals were modeled by a spherical infall model termed "The Great Attractor" (GA), 
with rig=A2 h'^Mpc toward (/, 6) = (309°, +18°), uzg = 535 kms'^ rc = 14.3 h'^Mpc, and 
n = 1.7 (Faber & Burstein 1988). A similar model vaguely flts the local infall of spirals into 
Virgo (Aaronson et al. 1982b) The merged Mark II data is well fltted by a multi-parameter 
hybrid consisting of a GA infall, a Virgo-centric infall, and a "local anomaly" - a bulk 
flow shared by the ~ 10 h~^Mpc-local neighborhood of 360 kms~^ perpendicular to the 
supergalactic plane. 

The toy models provide an intuitive picturing of the large-scale motions with clues 
about the associated mass sources, and they can be used as simple statistics for the 
comparison with theory (e.g. the bulk flow, §7.1; and the GA model, Bertschinger & 
Juszkiewicz 1988). However, toy modeling imposes an oversimplifled geometry associated 
with assumed "sources" on a complex velocity fleld which actually arises from a continuous 
fleld of asymmetric density fluctuations (§4.5, §5). Moreover, the bulk velocity statistic 
computed globally suffers from a bias due to the large-scale sampling gradients. The 
monopole, involving the radial decrease of sampling density and rise of errors, has the 
effect of reducing the effective volume and thus reducing the apparent conflict between the 
high bulk flow and the theoretical expectations (Kaiser 1988). The sampling dipole, arising 
for example from oversampling in the GA direction, tends to enhance the component of 
the bulk velocity in that direction because it is dominated by the velocity in a smaller 
effective volume (Regos & Szalay 1988; Szalay 1988). The sampling quadrupole arising 
from the Galactic zone of avoidance [zoa) introduces larger shot-noise into the component 
parallel to the Galactic plane, which could also result in an artiflcially high bulk velocity. 
These biases can could be partially cured by equal- volume weighting (§4.2). 

4.2. Potential Analysis 

If the LSS evolved according to GI, then the large-scale velocity fleld is expected 
to be irrotational, Vx« = 0. Any vorticity mode would have decayed during the linear 
regime as the universe expanded, and, based on Kelvin's circulation theorem, the flow 
remains vorticity-free in the quasi-linear regime as long as it is laminar. Irrotationality 
implies that the velocity fleld can be derived from a scalar potential, v{x) = — V$(s), so 
the radial velocity fleld u{x) should contain enough information for a full reconstruction. 
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In the POTENT procedure (Bertschinger & Dekel 1989) the potential is computed by 
integration along radial rays from the observer 

$(s) = - / u{r',e,(l))dr' , (19) 

Jo 

and the two missing transverse velocity components are then recovered by differentiation. 
Then 5[x) is approximated by 5c (eq. 7). The non-trivial step is the smoothing of the data 
into u[x). The aim in POTENT (Dekel et al. 1990) is to reproduce the u{x) that would 
have been obtained had the true v{x) been sampled densely and uniformly and smoothed 
with a spherical Gaussian window of radius Rg. With the data as available ui^x,.) is taken 
to be the value at x = Xc of an appropriate local velocity model v{Q,k,x — Xc) obtained by 
minimizing the sum (17) in terms of the parameters Q,k within an appropriate local window 
Wi = W{xi,Xc) chosen as follows. 

Tensor window. Unless Rg ^ r, the UiS cannot be averaged as scalars because the 
directions Xi differ from x^^ so u[Xc,) requires a fit of a local 3D model. The original 
POTENT used the simplest local model, v[x) =B, for which the solution can be expressed 
explicitly in terms of a tensor window function. 

Window bias. The tensorial correction to the spherical window has conical symmetry, 
weighting more heavily objects of large Xi-Xc. The resultant bias of a true infall transverse 
to the LOS is a flow towards the LG, e.g. ~ 300 kms~"^ at the GA in the current 
reconstruction. A way to reduce this bias is by generalizing B into a linear velocity model, 

v{x) = B + L- {x -Xc) , (20) 

with L a symmetric tensor, which ensures local irrotationality. The linear terms tend 
to "absorb" most of the bias, leaving v{Xc) = B less biased. Unfortunately, a high-order 
model tends to pick undesired small-scale noise. The optimal compromise is found to be a 
flrst-order model flt out to r = 40 h~"^Mpc, smoothly changing to a zeroth-order flt beyond 
60 h-^Mpc (Dekel et al. 1994). 

Sampling -gradient bias (SG). If the true velocity fleld is varying within the 
effective window, the non-uniform sampling introduces a bias because the smoothing is 
galaxy-weighted whereas the aim is equal-volume weighting. One should weight each 
object by the local volume it "occupies", e.g. Vi oc R^ where R^ is the distance to the 
n-th neighboring object {e.g. n = 4). This procedure is found via simulations to reduce 
the SG bias in Mark III to negligible levels typically out to 60 h~"^Mpc but away from the 
ZOA. The Rn{x) fleld can serve later as a flag for poorly sampled regions, to be excluded 
from any quantitative analysis. 

Reducing random errors. The ideal weighting for reducing the effect of Gaussian noise 
has weights Wi oc a~ but this spoils the carefully- designed volume weighting, biasing u 
towards its values at smaller r^ and at nearby clusters where the errors are small. A 
successful compromise is to weight by both, i.e. 

W{x,,Xc)^V^(j-^ exp[-{x,-Xcf/2Rl] . (21) 
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Note that POTENT could alternatively vary Rg to keep the random errors at a constant 
level, but at the expense of producing fields not directly comparable to theoretical models 
of uniform smoothing. Another way to reduce noise is to eliminate badly-observed galaxies 
with large residuals \ui — u[xi)\^ where u[xi) is obtained in a pre-POTENT smoothing stage. 
The whole smoothing procedure has been developed with the help of carefully- designed 
mock catalogs "observed" from simulations. 

Estimating the random errors. The errors in the recovered fields are assessed by 
Monte-Carlo simulations, where the input distances are perturbed at random using a 
Gaussian of standard deviation cr^ before being fed into POTENT. The error in (^ at a 
grid point is estimated by the standard deviation of the recovered 5 over the Monte-Carlo 
simulations, as (and similarly (Ty). In the well-sampled regions, which extend in Mark III 
out to 40—60 h~"^Mpc, the errors are as sa 0.1—0.3, but they may blow up in certain regions 
at large distances. To exclude noisy regions, any quantitative analysis should be limited 
to points where a^ and as are within certain limits. 

Several variants of POTENT are worth mentioning. The potential integration is 
naturally done along radial paths using only u[x) because the data are radial velocities, 
but recall that the smoothing procedure determines a 3D velocity using the finite 
effective opening angle of the window ~ Rs/r. The transverse components are normally 
determined with larger uncertainty but the non-uniform sampling may actually cause the 
minimum-error path to be non-radial, especially in regions where Rg/r^ 1. For example, 
it might be better to reach the far side of a void along its populated periphery rather than 
through its empty center. The optimal path can be determined by a max-flow algorithm 
(Simmons et al. 1994). It turns out in practice that only little can be gained by allowing 
non-radial paths because large empty regions usually occur at large distances where the 
transverse components are very noisy. Still, it is possible that the derived potential can be 
somewhat improved by averaging over many paths. 

The use of the opening angle to determine the transverse velocities can be carried one 
step further by fitting the data to a power-series generalizing the linear model, 

where X = x—Xc- If the matrices are all symmetric then the velocity model is automatically 
irrotational, and can therefore be used as the final result without appealing to the potential, 
and with the density being automatically approximated by 8(.{x) = \\I — Lij\\ — 1 (eq. 
7). The fit must be local because the model tends to blow up at large distances. The 
expansion can be truncated at any order, limited by the tendency of the high-order terms 
to pick up small-scale noise. The smoothing here is not a separate preceding step, the SG 
bias is reduced, and there is no need for numerical integration or differentiation, but the 
effective smoothing is again not straightforwardly related to theoretical models of uniform 
smoothing (Blumenthal & Dekel, in prep.). 

Another method of potential interest without a preliminary smoothing step is based 
on wavelet analysis, which enables a natural isolation of the structure on different scales 
(Rauzy et al. 1993). The effective smoothing involves no loss of information and no specific 
scale or shape for the wavelet, and the analysis is global and done in one step. How 
successful this method will be in dealing with noisy data and its comparison with theory 
and other data still remains to be seen. 
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4.3. Regularized Multi-Parameter Models — Wiener Filter 

A natural generalization of the toy models of §4.1 is a global expansion of the fields 
in a discrete set of basis functions, such as a Fourier series of the sort 

k 

8 = Y^[aj^sin(fc • x) + b^cos{k ■ x)], v = Y^ — [a^cos{k ■ x) — b^s'm{k ■ x)] , (23) 

with a very large number [2m) of free Fourier coefficients, and where 5 and v are related 
via the linear approximation (4). The maximization of the likelihood involves a 2m X 2m 
matrix inversion and, even if 2m is appropriately smaller than the number of data, the 
solution will most likely follow Murphy's law and blow up at large distances, yielding 
large spurious fluctuations where the data are sparse, noisy and weakly constraining, so 
the global flt requires some kind of regularization. Kaiser & Stebbins (1991) proposed 
to maximize the probability of the parameters subject to the data and an assumed prior 
model for the probability distribution of the Fourier coefficients; they assumed a Gaussian 
distribution with a PS (a? ) = (6? ) = Po^"- 

This is in fact an application of the Wiener filter method {e.g. Rybicki & Press 1992) 
for flnding the optimal estimator of a fleld 8{x) which is a linear functional of another fleld 
u(s), given noisy data Ui of the latter and an assumed prior model. It is 

8opt{x) = {8{x)ui) {u,Uj)~^ Uj , (24) 

where the indices run over the data (Hoffman 1994; Stebbins 1994). If u^ = u{Xi)-\-ei with e^ 
independent random errors of zero mean, then the cross-correlation terms are [8{x)u{xi)), 
and the auto-correlation matrix is 

{u^Uj) = {u{xi)u{xj)) + e]8,j , (25) 

both given by the prior model. 8opt is thus determined by the model where the errors 
dominate, and by the data where the errors are small. If the assumed prior is Gaussian, the 
optimal estimator is also the most probable fleld given the data, which is a generalization 
of the conditional mean given one constraint {e.g. Dekel 1981). This conditional mean 
fleld can be the basis for a general algorithm to produce constrained random realization 
(Hoffman and Ribak 1991). The same technique can also be applied to the inverse problem 
of recovering the velocity from observed density. Note (Lahav et al. 1994) that for Gaussian 
flelds the above procedure is closely related to maximum-entropy reconstruction of noisy 
pictures (Gull & Daniell 1978). 

The maximum-probability method has so far been applied in a preliminary way to 
heterogeneous data in a box of side 200 h~"^Mpc (to eliminate periodic boundary effects) 
with 18 h~"^Mpc resolution. No unique density fleld came out as different priors (—3 < n < 1) 
led to different flts with similar x^ ■ This means that the data used were not of sufflcient 
quality to determine the flelds this way and it would be interesting to see how this method 
does when applied to better data. The method still lacks a complete error analysis, it 
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needs to somehow deal with non-hnear effects, and it needs to correct for IM bias (perhaps 
as in §3.2). 

A general undesired feature of maximum probability solutions is that they tend to be 
oversmoothed in regions of poor data, relaxing to (^ = in the extreme case of no data. 
This is unfortunate because the signal of true density is modulated by the density and 
quality of sampling - a sampling bias which replaces the SG bias. The effectively varying 
smoothing length can affect any dynamical use of the reconstructed field {e.g. deriving 
a gravitational acceleration from a density field), as well as prevent a straightforward 
comparison with other data or with uniformly-smoothed theoretical fields. Yahil (1994) 
has recently proposed a modified filtering method which partly cures this problem by 
forcing the recovered field to have a constant variance. 

4.4. Malmquist-Free Analysis 

The selection bias (§3.2) can be practically eliminated from the calibration of an 
inverse TF relation, rji^M), as long as the internal velocity parameter rj does not explicitly 
enter the selection process. An inverse analysis requires assuming a parametric model for 
the velocity field, v^a^k^x) (Schechter 1980). The sum minimized instead of (17) is 



/ Wi [rii{ohserved) — rji (model)] , (26) 



with the model r] given by the inverse TF relation, 

rii[model) = a + bMi = a + b[mi — Slogr^) , ri = Zi — Xi ■ w(ayt, Xi) . (27) 

The parameters are the inverse TF parameters a and b and the a,t characterizing the 
velocity model. This method is indeed free of inferred-distance M bias as long as r^ is 
uniquely derived from Zi^ and not from the inferred diatance. 

An inverse method was first used by Aaronson et al. (1982b) to fit a Virgo-centric 
toy model to a local sample of spirals, and attempts to implement the inverse method 
to an extended sample with a general velocity model are in progress (MFPOT by Yahil 
et al. 1994). This is a non-trivial problem of non-linear multi-parameter minimization 
with several possible routes. If the velocity model is expressed in 2;-space, then r^ is given 
explicitly by Zi but the results suffer from oversmoothing in collapsing regions. If in r-space, 
then ri is implicit in the second equation of (27) requiring iterative minimization, e.g. by 
carrying r^ from one iteration to the next one. The velocity model could be either global 
or local, with the former enabling a simultaneous minimization of the TF and velocity 
parameters and the latter requiring a sequential minimization of global TF parameters 
and local velocity parameters. While the results are supposed to be free of Malmquist 
bias, they suffer from other biases which have to be carefully diagnosed and corrected for. 
The inverse results will have to The inverse method is also being generalized to account for 
small-scale velocity noise (Willick, Burstein & Tornen 1994). The Malmquist-free results 
will have to be consistent with the M-corrected forward results before one can put to rest 
the crucial issue of M bias and its effect on our LSS results. 
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4.5. Fields of Velocity and Mass Density 

Figure 3 shows supergalactic-plane maps of the velocity field in the CMB frame and the 
associated 8c field (for O = 1) as recovered by POTENT from the preliminary Mark III data. 
The data are reliable out to ~ 60 h~"^Mpc in most directions outside the Galactic plane 
(y = 0), and out to ~ 70 h~"^Mpc in the direction of the GA (left-top) and Perseus-Pisces 
(PP, right-bottom). Both large-scale (~ 100 h~"^Mpc) and small-scale (~ 10 h~"^Mpc) 
features are important; e.g. the bulk velocity refiects properties of the initial fiuctuations 
and of the DM (§7), while the small-scale variations indicate the value of O (§8). 

The velocity map shows a clear tendency for motion from right to left, in the general 
direction of the LG-CMB motion (L, 5 = 139°, -31°). The bulk velocity within 60 h'^Mpc 
is 300 — 350 km s,~^ towards (L, B sa 166°, —20°) (§7.1) but the fiow is not coherent over the 
whole volume sampled, e.g. there are regions in front of PP and at the back of the GA where 
the XY velocity components vanish, i.e. the streaming relative to the LG is opposite to 
the bulk fiow direction. The velocity field shows local convergences and divergences which 
indicate strong density variations on scales about twice as large as the smoothing scale. 

The GA at 12 h~"^Mpc smoothing and O = 1 is a broad density peak of maximum height 
(^=1.2 lb 0.3 located near the Galactic plane y = at X sa — 40 h~"^Mpc. The GA extends 
towards Virgo near V sa 10 (the "Local Supercluster"), towards Pavo-Indus-Telescopium 
(PIT) across the Galactic plane to the south {Y < 0), and towards the Shapley 
concentration behind the GA. The structure at the top roughly coincides with the "Great 
Wall" of Coma, with 8 sa 0.5. The PP peak which dominates the right-bottom is centered 
near Perseus with (^= 1.0 it 0.4. PP extends towards Aquarius in the southern hemisphere, 
and connects to Cetus near the south Galactic pole, where the "Southern Wall" is seen in 
redshift surveys. Underdense regions separate the GA and PP, extending from bottom-left 
to top-right. The deepest region in the Supergalactic plane, with 8 = —0.7 it 0.2, roughly 
coincides with the galaxy-void of Sculptor (Kauffman et al. 1991). 

One can still find in the literature statements questioning the very existence of the 
GA {e.g. Rowan- Robinson 1993), which simply refiect ambiguous definitions for this 
phenomenon. A GA clearly exists in the sense that the dominant feature in the local 
inferred velocity field is a coherent convergence, centered near X sa —40. It is another 
question whether the associated density peak has a counterpart in the galaxy distribution 
or is a separate, unseen entity. The GA is ambiguous only in the sense that the good 
correlation observed between the mass density inferred from the velocities and the galaxy 
density in redshift surveys is perhaps not perfect (§6.2). 

Other cosmographic issues of debate are whether there exists a back-fiow behind the 
GA in the CMB frame, and whether PP and the LG are approaching each other. These 
effects are detected by the current POTENT analysis only at the 1.5(T level in terms of the 
random uncertainty. Furthermore, the freedom in the zero-point of the distance indicators 
permits adding a Hubble-like peculiar velocity which can balance the GA back fiow and 
make PP move away from the LG. Thus, these issues remain debatable. 

To what extent should one believe the recovery in the ZOA which is empty of tracers? 
The velocities observed on the two sides of the ZOA are used as probes of the mass in 
the ZOA. The interpolation is based on the assumed irrotationality, where the recovered 





Figure 3: The fluctuation flelds of velocity and rriass- density in the Supergalactic 
plane as recovered by POTENT from the Mark III velocities of ~ 3000 galaxies 
with 12 h~"^Mpc smoothing. The vectors shown are projections of the 3D velocity 
fleld in the CMB frame. Distances and velocities are in 1000 kms~"^. Contour 
spacing is 0.2 in 8, with the heavy contour marking 8 = Q and dashed contours 
8 < Q. The LG is at the center. The GA is on the left, PP on the right, and Coma 
is at the top. The grey-scale in the contour map and the height of the surface in 
the landscape map are proportional to 8 (Dekel et al. 1994). 
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transverse components enable a reconstruction of the mass-density. However, while the SG 
bias can be corrected where the width of the ZOA is smaller than the smoothing length, 
the result could be severely biased where the unsampled region is larger. With 12 h~"^Mpc 
smoothing in Mark III the interpolation is suspected of being severely biased in ~ 50% of 
the ZOA at r = 40 h~"^Mpc (where i?4 > -R^), but the interpolation is pretty safe in the 
highly populated GA region, for example. Indeed, a deep survey of optical galaxies at small 
Galactic latitudes recently discovered a very big cluster centered at (/, 6, z) = (325°, —7°, ~ 
45 h~"^Mpc) (R. Kraan-Korteweg, private comm.), very near the central peak of the GA 
as predicted by POTENT, ~(320°,0°,40 h'^Mpc) (Kolatt et al. 1994). 

5. PREDICTED MOTIONS FROM THE GALAXY ^-DISTRIBUTION 

All-sky complete redshift surveys provide extremely valuable data complementary 
to the peculiar velocity data, and efficient techniques for measuring redshifts make 
them deeper, denser, and more uniform. Under the assumption of GI and an assumed 
biasing relation between galaxies and mass, a redshift survey enables an independent 
reconstruction of the density and velocity fields. The comparison of the smoothed fields 
recovered from redshifts and from velocities is a very important tool for testing the basic 
hypotheses and for determining the cosmological parameters. 

The solution to the linearized GI equation V-t> = —f8 for an irrotational field is 

v{x) = — f d'x'5{x') ^[2''i3 ■ (28) 

^'^ ./all space P — X\ 

The velocity is proportional to the gravitational acceleration, which ideally requires full 
knowledge of the distribution of mass in space. In practice (Yahil et al. 1991) one is 
provided with a flux-limited, discrete redshift survey, obeying some radial selection function 
(j){r). The galaxy density is estimated by 

1 + Sg{x) = J2 ^"V(r^)-^<^Lac(^ - ^^) , (29) 

where n = V~^ ^(l){ri)~^ is the mean galaxy density, and the inverse weighting by (/> 
restores the equal-volume weighting. Eq. (28) is then replaced by 



v[x] = —l d^x' 8Jx')S{\x' -x\)^, % . (30) 



47r -/r<R 



Under the assumption of linear biasing, 8g = h8, the cosmological dependence enters through 
fi = f[Q)/b. The integration is limited to r < Rmax where the signal dominates over 
shot-noise. S{'y) is a small-scale smoothing window (>500 kms~"^) essential for reducing 
the effects of non-linear gravity, shot-noise, distance uncertainty, and triple-value zones. 
The distances are estimated from the redshifts in the LG frame by 

r, = z,-x,-[v{Xi)-v{{))] . (31) 

Equations (30-31) can be solved iteratively: make a flrst guess for the s^, compute the Vi 
by eq. (30), correct the Xi by eq. (31), and so on until convergence. The convergence can 
be improved by increasing fi gradually during the iterations. Relevant issues follow. 
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Selection function. An accurate knowledge of the probability that a galaxy at a given 
distance be included in the sample is essential, especially at large distances where c/)""^ 
can introduce large errors. For a given flux limit, (j) can be evaluated together with 
the luminosity function using a maximum-likelihood technique independent of density 
inhomogeneities . 

Zone of Avoidance. Regions in the sky not covered by the survey have to be flUed 
with mock galaxies by some method of extrapolation from nearby regions. One way is to 
distribute these galaxies Poissonianly with the mean density of an adjacent volume, or to 
actually clone the adjacent region (Hudson 1993a). A more sophisticated extrapolation 
uses spherical harmonics and the Wiener fllter method {e.g. Lahav et al. 1994). 

Triple-valued zones. Galaxies in three different positions along a LOS through a 
contracting region may have the same redshift. Given a redshift in a collapsing region 
where the problem is not resolved by the smoothing used, one can either take some average 
of the three solutions, or make an intelligent choice between them, e.g. by using the velocity 
fleld derived from observed velocities. 

Estimating shot-noise. This major source of error due to the flnite number of galaxies 
can be crudely estimated using bootstrap simulations, where each galaxy is replaced with k 
galaxies, k being a Poisson deviate of (k) = 1. For each realization one calculates (j) and n, 
corrects for the ZOA^ and solves for the linear velocity fleld. The mean and variance of the 
resulting density fleld are measures of the systematic and random errors. The bootstrap 
simulations demonstrate that the uncertainty in 5g from the 1.2 Jy IRAS sample is typically 
less than 50% of the uncertainty in the density derived by POTENT from observed Mark 
III velocities. 

Nonlinear Biasing. Galaxies need not be faithful tracers of the mass {e.g. Dekel 
& Rees 1987), but there is growing evidence that they are strongly correlated (§6.2). 
This correlation can be crudely assumed to be a deterministic relation between the local 
smoothed density flelds, e.g. linear biasing 5g = bS^ which is one realization of the linear 
statistical relation between the variances of the flelds predicted for linear density peaks in a 
Gaussian fleld (Kaiser 1984; Bardeen et al. 1986). However, a more sophisticated analysis 
may require a more realistic biasing relation, e.g. an exception from linear biasing which 
must be made for negative 5g and 6 < 1 to prevent 5 from falling unphysically below —1. 
The non-linear generalization 1 -\- Sg = {l -\- S) is useful {e.g. Dekel et al. 1993) and it flts 
quite well the biasing seen in simulations of galaxy formation in a CDM scenario involving 
cooling and gas dynamics (Cen & Ostriker 1993), except that a small correction is needed 
to force the means of 5 and 5g to vanish simultaneously as required by deflnition. 

Quasilinear correction. Even after 12 h~"^Mpc smoothing, 5g is of order unity in places, 
necessitating a quasi-linear treatment. Local approximations from v to 5 were discussed 
in §2, but the non-local nature of the inverse problem makes it less straightforward. A 
possible solution is to flnd an inverse relation of the sort V'i> = -F(0, 5g)^ including non-linear 
biasing and non-linear gravity. This is a Poisson-like equation in which —fi8g{x) is replaced 
by -F(x), and since the smoothed velocity fleld is irrotational for quasi-linear perturbations 
as well, it can be integrated analogously to eq. (30). With smoothing of 10 h~"^Mpc and 
(3 = 1, the approximation based on 5^ has an rms error < 50 kms~"^ (Mancinelli et al. 
1994). Note that for very small b the 5 associated with the observed 5g could be non-linear 
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Figure 4: The fluctuation flelds of galaxies in the Supergalactic plane as deduced 
from redshift surveys with 10 h~"^Mpc smoothing. Distances and velocities are 
in kms~"^. Density contour spacing is 0.2. (a) Reconstructed by A. Yahil and 
M. Strauss using a power-preserving fllter from the IRAS 1.2 Jy data in a sphere 
of radius 160 h~"^Mpc. (b) Reconstructed by Hudson (1993a, 1994) from optical 
data within 80 h~"^Mpc, extrapolated into the unsampled areas ouside the heavy 
contour. r^r^ 



to the extent that the quasi-hnear approximations break down. When 8 is given, the 
Poisson equation can be integrated more efficiently by using grid-based FFT techniques 
than by straightforward summation. The r-space 8{x) deduced from the 2;-space galaxy 
distribution is not too sensitive to non-linear effects or the value of O, so it is a reasonable 
shortcut to correct for non-linear effects only in the transformation from 8 to v using FFT. 

The IRAS Point Source Catalog served as the source for two very valuable redshift 
surveys, which have been carried out and analyzed in parallel. One contains the 5313 
galaxies brighter than 1.2 Jy at 60/i with sky coverage (almost) complete for |6| > 5° 
covering 88% of the sky (Strauss et al. 1990a, 1992, to 1.9 Jy; extended by Fisher 1992). 
The other is a l-in-6 sparsely-sampled survey of sa 2300 galaxies down to the IRAS 
flux limit of 0.6 Jy (Rowan- Robinson et al. 1990, QDOT), which is now being extended 
to a fully-sampled survey (Saunders et al. in prep.). As for optical galaxies, Hudson 
(1993a) developed a clever way to reconstruct a statistically uniform density fleld out to 
~ 80 h~"^Mpc by combining the UGC/ESO diameter-limited angular catalogs and the 
ZCAT incomplete redshift survey. Figure 4 shows maps of the galaxy density flelds and 
the associated predicted velocity flelds, with the main features corresponding to those 
recovered from observed velocities (Fig. 3), e.g. the GA, PP and Coma superclusters and 
the voids in between (§6.2). 

6. TESTING BASIC HYPOTHESES 

The data and analyses described above, combined with other astrophysical data, can 
help us evaluate some of the basic hypotheses laid out in the Introduction. 

6.1. CMB Fluctuations versus Motions 

If the CMB deflnes a standard cosmological frame then the established helio-centric 
dipole pattern of 8T jT = (1.23ib0.01) X 10~^ is a direct measurement of peculiar velocity of 
theLG: 1^(0) = 627 ± 22 km s'^ towards (/, 6) = (276° ± 3°, +30° ± 3°) (Kogut et al. 1993). 
The Copernican hypothesis then implies that large peculiar velocities exist in general, and 
the question is only how coherent are they. A more esoteric interpretation tried to explain 
the CMB dipole by a horizon-scale gradient in entropy, relic of bubbly inflation, which 
can be made consistent with the smallness of the quadrupole and the achromaticity of the 
dipole {e.g. Gunn 1988; Paczynski & Piran 1990). In addition to the general objection 
to global anisotropy based on the simplicity principle of Occam's razor, a non-velocity 
interpretation fails to explain the fact that the gravitational acceleration vector at the LG, 
5'(0), as inferred from the galaxy distribution in our cosmological neighborhood (§8.1), is 
within 20° of V(0) and of a similar amplitude of several hundred km s,~^ for any reasonable 
choice of O. This argument in favor of the reality of V(0) and the standard GI picture 
is strengthened by the fact that a similar g{Q) is obtained from the POTENT mass fleld 
derived from velocity divergences in the neighborhood (Kolatt et al. 1994). 

The measurements of CMB fluctuations on scales < 90° are independent of the 
local streaming motions, but GI predicts an intimate relation between them. The CMB 
fluctuations are associated with fluctuations in gravitational potential, velocity and density 
in the surface of last scattering at 2; ~ 10^, while similar fluctuations in our neighborhood 
have grown by gravity to produce the dynamical structure observed. The comparison 
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between the two is therefore a crucial test for GI. Before COBE, the streaming velocities 
were used to predict the expected level of CMB fluctuations. The local surveyed region of 
~ 100 h~"^Mpc corresponds to a ~ 1° patch on the last-scattering surface. The major effect 
on scales > 1° is the Sachs- Wolfe effect (1967), where potential fluctuations A$(, induce 
temperature fluctuations via gravitational redshift, 

5T A$g Vx . , 

Since the velocity potential is proportional to $(, in the linear and quasi-linear regimes, 
A$(, is ~ Vx^ where x is the scale over which the bulk velocity is V. A typical bulk 
velocity of ~ 300 kms~"^ across ~ 100 h~"^Mpc corresponds to 5T/T ^ 10~^ at ~ 1°. If 
the fluctuations are scale-invariant (n = 1), then ST/T ~ 10~^ is expected on all scales 
> 1°. Bertschinger et al. (1990) produced a 5T/T map of the local region as seen by a 
distant observer, and predicted 5T/T ^ 10~^ from the local potential well associated with 
the GA. Now that CMB fluctuations of sa 10~^ have been detected in the range 1° — 90° 
{e.g. Smoot et al. 1992; Schuster et al. 1993; Cheng et al. 1993), the argument can be 
reversed: if one assumes GI then the expected bulk velocity in the surveyed volume is 
~300 kms~"^, i.e. the motions are real. If, alternatively, one accepts the velocities as real, 
then the CMB-POTENT agreement is relatively a sensitive test of GI, truly addressing 
the speciflc time- evolution of structure predicted by GI. 

6.2. Galaxies versus Dynamical Mass 

The theory of GI plus the assumption of linear galaxy biasing predicts a correlation 
between the dynamical density fleld and the galaxy density fleld, which can be addressed 
quantitatively based on the estimated errors in the two data sets. Figure 5 compares density 
maps in the Supergalactic plane for IRAS 1.2 Jy galaxies and POTENT Mark III mass, 
both Gaussian smoothed with radius 12 h~"^Mpc. The correlation is evident - the GA, PP, 
Coma and the voids all exist both as dynamical entities and as structures of galaxies. A 
quantitative comparison of these new data is in progress, but so far an elaborate statistical 
analysis (Dekel et al. 1993) has been applied only to the earlier POTENT reconstruction 
based on Mark II data and IRAS 1.9 Jy survey. Noise considerations in POTENT limited 
that analysis to a volume ~ (53 h~"^Mpc)^ containing ~ 12 independent density samples. 
Monte-Carlo noise simulations showed that the POTENT density is consistent with being 
a noisy version of the IRAS density, i.e. the data are in agreement with the hypotheses of 
GI plus linear biasing. 

What is it exactly that one can learn from the observed correlation (Babul et al. 
1994)? First, it is hard to invoke any reasonable way to make the galaxy distribution and 
the TF measurements agree so well unless the velocities are real (HG), provided that IM 
bias has been properly corrected (§3.2, §4.4). Then, it is true that gravity is the only 
long-range force which could attract galaxies to stream toward density concentrations, but 
the fact that this correlation is predicted by GI plus linear biasing does not necessarily 
mean that it can serve as a sensitive test for either. Recall that converging (or diverging) 
flows tend to generate density hills (or valleys) simply as a result of mass conservation, 
independent of the source of the motions. 
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Figure 5: POTENT mass versus IRAS galaxy density fields in the Supergalactic 
plane, Gaussian smoothed with radius 12 h~"^Mpc. Contour spacing is 0.2 (Dekel 
et al. 1993). 

Let us assume for a moment that galaxies trace mass, i.e. the linearized continuity 
equation, 8 = —V • «, is valid for the galaxies as well. The observed correlation means 
(^oc— V-w, and together they imply that 8(x8, or equivalently that V-w is proportional to 
its time average. This property is not exclusive to GI; one can construct a counterexample 
where the velocities are produced by a non-GI impulse. Even irrotationality does not 
follow from 8 oc — V-t>; it has to be adopted based on theoretical arguments in order to 
enable reconstruction from radial velocities or from observed densities. Once continuity 
and irrotationality are assumed, the observed 8(x — V-t> implies a system of equations which 
is identical in all its spatial properties to the equations of GI, but can differ in the constants 
of proportionality and their temporal behavior. It is therefore impossible to distinguish 
between GI and a non-GI model which obeys continuity plus irrotationality based only on 
snapshots of present-day linear fluctuation flelds. This makes the relation between CMB 
fluctuations and velocities an especially important test for GI. On the other hand, the 
fact that the constant of proportionality in (^ oc — V-w is indeed the same everywhere is a 
non-trivial requirement from a non-GI model. A version of the explosion scenario (Ostriker 
& Cowie 1981; Ikeuchi 1981), which tested successfully both for irrotationality and v — 8 
correlation, requires certain synchronization among the explosions (Babul et al. 1994). 

While the sensitivity of the v — 8 relation to GI is only partial, this relation turns out 
to be quite sensitive to the validity of a continuity-like relation for the galaxies; when the 
latter is strongly violated all bets are off for the 'V—8g relation. A non-linear biasing scheme 
would make continuity invalid for the galaxies, which would ruin the t}—8g relation even if 
GI is valid. The observed correlation is thus a sensitive test for density biasing. It implies, 
subject to the errors, that the ~ 10 h~"^Mpc-smoothed densities of galaxies and mass are 
related via a simple, nearly linear biasing relation with b not far from unity. 
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6.3. Environmental Effects: Ellipticals versus Spirals 

A priori, it is possible that the inferred motions are just a reflection of systematic 
variations in the distance indicators, e.g. due to environmental variations in intrinsic 
galaxy properties, or in the apparent quantities due to Galactic absorption {e.g. Silk 1989; 
Djorgovski et al. 1989). Efforts to detect correlations between velocities and certain galaxy 
or environmental properties have led so far to null or at most minor detections (Aaronson 
& Mould 1983; Lynden-BeU et al. 1988; Burstein 1990a, 1990b; Burstein et al. 1990) - 
no correlation with absorption and with local galaxy density, and marginal correlations 
with absolute luminosity and with stellar population (Gregg 1993). Admittedly, these null 
results are only indicative as one cannot rule out a correlation with some other property 
not yet tested for. 

Qualitative comparisons of the velocities of ellipticals (Es) and spirals (Ss) indicated 
general agreement (Burstein et al. 1990), and now the Mark III data enable a quantitative 
comparison at the same positions (Kolatt & Dekel 1994a). Figure 6 compares the two 
flelds as derived independently by POTENT, showing a general resemblance. The radial 
velocities of each type were interpolated into a smoothed fleld u{x) on a grid using 
12 h~"^Mpc POTENT smoothing (§4.2), and then compared within a volume limited by the 
poor sampling of Es to ~ (50 h~"^Mpc)^, containing ~ 10 independent sub-volumes. The 
two flelds were found to be consistent with being noisy versions of the same underlying fleld, 
while the opposite hypothesis of complete independence is strongly ruled out, indicating 
that the consistency is not dominated by errors. A possible discrepancy indicated earlier 
by the Mark II data (Bertschinger 1991) is now gone with the improved data and 
bias-correction. The strength of this test will improve further as the E sample grows in 
extent {e.g. project EFAR by CoUess et al. 1993) and in accuracy {e.g. the SBF method). 
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Figure 6: Ellipticals versus Spirals. Mass-density and velocity flelds in the 
Supergalactic plane as recovered by POTENT independently for the two types. 
Smoothing is 12 h~"^Mpc. Contour spacing is 0.2. Distances and velocities are in 
1000 kms-i (Kolatt & Dekel 1994). 
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The E-S correlation is consistent with motions (H6), but it cannot rule out 
environmental effects beyond making this idea less plausible. If such effects were dominant, 
then the two distance indicators would have to vary coherently, which would require 
that the large-scale properties of the environment be different from the local properties 
determining the galaxy type (E in clusters, S in field). It would also require that the 
properties affecting the virial equilibrium which determines the Dn—(T relation, and those 
affecting the additional constraint involved in the TF relation at galaxy formation (§2.1), 
varied together in space. The one quantity that could plausibly affect the two in a 
correlated way is M/L, but to test for this one will have to compare the TF/D„— cr 
velocities to those inferred from an independent indicator, e.g. SBF. 

The E-S correlation is consistent with the hypothesis of GI (114) because, following 
Galileo, the velocities of all test bodies in a given gravitational potential are predicted to 
be the same as long as they share the same initial conditions. However, the observed E-S 
correlation does not rule out any non-GI model where all objects obtain the same velocities 
independent of their type, e.g. cosmological explosions or radiation pressure instabilities 
(Hogan & White 1986). 

A practical use of the independent derivation of u[x) for Ss and Es is in matching the 
zero points of the distance indicators, otherwise determined in an arbitrary way (§2). A 
5% Hubble-like outflow has to be added to the peculiar velocities of Es from Mark II to 
optimally match the S data of Mark III. 

7. THE INITIAL FLUCTUATIONS 

Having assumed GI, the LSS can be traced backward in time to recover the initial 
fluctuations and constrain statistics which characterize them as a random fleld, e.g. the 
power spectrum (PS), and the probability distribution functions (PDFs). "Initial" here 
may refer either to the linear regime at z ^ 10^ after the onset of the self-gravitating matter 
era, or to the origin of fluctuations in the early universe before being flltered on sub-horizon 
scales during the plasma-radiation era. The PS is flltered on scales < 100 h~"^Mpc by 
DM-dominated processes, but its shape on scales > 10 h~"^Mpc is not affected much 
by recent non-linear effects (because the rapid density evolution in superclusters roughly 
balances the slow evolution in voids at the same wavelength). The shape of the one-point 
PDF, on the other hand, is expected to survive the plasma era unchanged but it develops 
strong skewness even in the mildly-non-linear regime. Thus, the present day PS can be 
used as is to constrain the origin of fluctuations (large scale) and the nature of the DM 
(small scale), while the PDF needs to be traced back to the linear regime flrst. (Review: 
Efstathiou 1990.) 

7.1 Power Spectrum — Dark Matter 

The competing LSS formation scenarios are reviewed in Peebles (1993, §25). If the 
DM is all baryonic, then by nucleosynthesis constraints (see Kolb & Turner 1990, §4) the 
universe must be of low density, O < 0.2, and a viable model for LSS is the Primordial 
Isocurvature Baryonic model (PIB) with several free parameters, typically of large relative 
power on large scales. With O ~ 1 the non-baryonic DM constituents are either "hot" or 
"cold", and the main competing models are CDM, HDM, and MDM - a 3:7 mixture of 
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the two {e.g. Blumenthal et al. 1988; Davis et al. 1992; Klypin et al. 1993). The main 
difference in the DM effect on the PS arises from free-streaming damping of the "hot" 
component of fluctuations on galactic scales. 

Bulk velocity. A simple and robust statistic related to the PS is the amplitude V of 
the vector average of the smoothed velocity fleld, «, over a volume deflned by a normalized 
window function WR{r) of a characteristic scale R {e.g. top-hat), 

(fxWR{x)v{x) , {V^) = ^ / dkP{k)Wl{k) . (33) 

{V^) is predicted for a linear model with a density spectrum P{k), where W'^{k), the 
Fourier transform of WR{r)^ emphasizes waves > R. The bulk velocity is obtained from 
the observed radial velocities by minimizing eq. (17). The report by Dressier et al. (1987) 
for the 7S Es sampled within ~60 h'^Mpc was 1^ = 599 ±104 towards (/, 6) = (312°, -^6°), 
which was interpreted prematurely as being in severe excess of the predictions of common 
theories. However, this measurement cannot be directly compared to the predictions for 
a top-hat sphere because the effective window is much smaller due to the nonuniform 
sampling (SG) and weighting (Kaiser 1988). The SG bias can be crudely corrected by 
volume weighting as in POTENT (§4.2), at the expense of large noise. Courteau et al. 
(1993) flnd for the tentative Mark III data: V40 = 466±48 (312°, +15°) and Veo = 398±47 
(305°, -|- 14°), where Vr refers to a top-hat sphere of radius R h~"^Mpc. Alternatively, V 
can be computed from the POTENT v fleld by simple vector averaging from the grid. 
Figure 7 shows two results, one minimizing the SG bias by Vi weighting, and the other 
reducing the random errors by weighting oc a~ (§4.2). Veo is found to be in the range 
270 — 360 kms~"^ (296°, -|-11°) - smaller than previous estimates. The additional random 
error from Monte-Carlo noise simulations is typically 15%, not including cosmic scatter 
due to the fact that only one sphere has been sampled. 

Mach number. The bulk velocity is robust but it relates to the normalization of the PS, 
not predicted from flrst principles by any of the competing theories but rather normalized 
by some other uncertain observation, e.g. the CMB fluctuations or the galaxy distribution 
with an unknown biasing factor. A statistic which measures the shape of the PS free of its 
normalization is the cosmic Mach number (Ostriker & Suto 1989), deflned as A^ = V/ S^ 
where S is the rms deviation of the local velocity from the bulk velocity, 

/x2 i-oo 
d'xWR{x)[v{x)-V]\ {S^) = ^ dkP{k)[l-Wl{k)]. (34) 
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Ai{R^ Rs) measures the ratio of power on large scales ^R to power on small scales ^ R 
Strauss et al. (1993) derived A^ = 1.0 for the local Ss (Aaronson et al. 1982a) with i? ~ 
20 h~"^Mpc and Rg — >0, and found ~5% of their CDM simulations to have Ai as large - 
a marginal rejection or consistency depending on taste. On larger scales, using POTENT 
with R = 60 h~"^Mpc (top-hat) and Rg = 12 h~"^Mpc (Gaussian), the tentative Mark III 
data yield A^ ~ 1, which roughly coincides with the rms expected from CDM but has 
only ~ 5% probability to be that low for a PIB spectrum (O = 0.1, h = 1, fully ionized, 
normalized to crs = l) (Kolatt et al. , in prep.). 
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Figure 7: The bulk velocity in a top-hat sphere of radius Reff about the 
LG as recovered by POTENT from the Mark III data. Shown are \V\ (filled), 
Vx (triangles), Vy (squares), and V^ (hexagons). Distances and velocities are 
in kms~"^. The two results shown reflect the systematic uncertainty. The la 
uncertainty due to random distance errors is ~ 15% (Dekel et al. 1994). 

Power spectrum. The velocity fleld by POTENT enables a preliminary determination 
of the mass PS itself in the range 10-100 h'^Mpc (Kolatt & Dekel 1994b). It is best 
determined by the potential fleld, which is smoother and less sensitive to non-linear effects 
than its derivatives v and 5. The PS was compared with the predictions of theoretical 
models which were N-body simulated, "observed", and fed into POTENT before the PS 
was computed using the same procedure applied to the observed data, thus eliminating 
systematic errors and also estimating random errors. The preliminary results indicate that 
the shape of the PS in the limited range sampled resembles a CDM-like PS with a shape 
parameter F ~ 0.5 (referring in CDM to 0/z), with the power index bending toward n ~ by 
100 h~"^Mpc. The normalization for the mass PS, also obtained by Seljak & Bertschinger 
(1994), is fas = 1.3 it 0.3. This is not too sensitive to the PS shape, but still note that it 
mostly determined by data at ~30 h~"^Mpc. With ergo ~ 1 for optical galaxies, this implies 
/3o ~ 1-3 lb 0.3 (compare §8). If O = 1, the quadrupole in 5T/T by COBE corresponds to 
ag = 1.0,0.6,0.5 with ~ 45% error (including cosmic scatter) for CDM, MDM and HDM 
spectra (yet to be calculated for an open universe), i.e. CDM is flne while HDM and MDM 
are ~2(T low. In comparison, the PS of the different luminous objects are all well described 
by a CDM-like PS with F sa 0.25, with the relative bias factors for Abell clusters, optical 
galaxies and IRAS galaxies in the ratios 4.5 : 1.3 : 1 (±6% rnis , e.g. Peacock & Dodds 
1994). If indeed aso ~ 1 for optical galaxies, then IRAS galaxies are slightly anti-biased 
and /3jSa 1.3/3o- 

A discrepancy on a very large scale? The galaxy velocity fleld in the local ~60 h~"^Mpc 
which has been studied quite accurately seems to be in general agreement with other 
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observations and with the predictions of the common GI scenarios, but there is a 
disturbing hint for a possible discrepancy with our basic hypotheses on larger scales. Lauer 
and Postman (1993, LP) measured the bulk flow of the volume-limited system of 119 
Abell/ACO clusters within z < 150 h~"^Mpc, estimating distances from luminosity versus 
slope of the surface-brightness proflle in brightest cluster galaxies (BCG) with claimed 16% 
error. The LG motion with respect to this system is found to be 561 it 284 kms~"^ towards 
(/,6) = (220°,-28°), roughly 80° off the direction of the LG-CMB velocity (276°, +30°). 
The inferred bulk velocity of the cluster system (of effective radius ~ 80 — 110 h~"^Mpc) 
relative to the CMB is F = 689 ± 178 kms'^ towards (343°, -^52°) (±23°), in the general 
direction of the LG motion, the galaxy bulk flow, the GA and the background Shapley 
concentration. Taken at face value this is a big velocity over a big scale; the rms bulk 
velocity within R= 100 h~"^Mpc as predicted by the common theories normalized to COBE 
is below 200 kms~"^. However, the errors are big too! 

The data have been subjected by LP to careful statistical tests. With a few tens of 
points contributing to the velocity in each direction, the shot-noise can clearly contribute a 
false signal of several hundred km s~^ . For assessing the theoretical implications, the actual 
observing scheme of LP was applied to N-body simulations of several competing scenarios 
(Strauss et al. 1994). Clusters were placed at the 119 highest density peaks of appropriate 
mass, and Galactic extinction and observational errors were modeled. As noted by LP, 
the fact that the measured velocity vector lies away from the ZOA increases the statistical 
signiflcance of this measurement. Taking the error ellipsoid into account, the probability 
by means of x^ of the LP result in the model universes simulated is 2.6-5.8%. Feldman 
& Watkins (1994), using a different technique, flnd the probabilities to be 6-10%. The 
LP bulk velocity is thus a ~2(T deviation from several common theories. This is probably 
not enough for a serious falsiflcation of these models, but it is certainly an intriguing 
result which motivates more accurate investigations. The BCG method is limited to one 
measurement per cluster, so an obvious strategy to reduce the errors would be to collect 
many TF and Dn—(T distances per cluster. 

A true ~ 700 kms~"^ velocity at i?~ 100 h~"^Mpc would be in serious conflict with GL 
First, if this velocity is typical then it predicts larger CMB fluctuations than observed on 
~2°. Second, the gravitational acceleration on the LP sphere as estimated from the spatial 
distribution of clusters on even larger scales (§8.1, Scaramella et al. 1991) predicts a flow 
of only ~200 kms~"^. The way to interpret the LP result in the context of the conventional 
theories is either as a ~ 2— 2.5(T statistical fluke, or as a biased result due to a yet-unresolved 
systematic error in the BCG method or in the sample, e.g. a signiflcant velocity biasing of 
the clusters (which is contrary to the naive expectations from GI, however). 

7.2 Back in Time 

The forward integration of the GI equations by analytic approximations or by N-body 
simulations cannot be simply reversed despite the time reversiblity of gravity. It is 
especially hopeless in collapsed systems where memory has been erased, but the case is 
problematic even for linear systems. When attempting backwards integration, the decaying 
modes {e.g. Peebles 1993, §5), having left no detectable trace at present, to, would amplify 
noise into dominant spurious fluctuations at early times. This procedure has a negligible 
probability of recovering the very special initial state of almost uniform density and tiny 
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velocities which we assume for the real universe at ti. The problem is of mixed boundary 
conditions: some of the six phase-space variables per particle are given at ti and some at 
to. This problem can be solved either by eliminating the decaying modes or by applying 
the principle of least action. 

Zel'dovich time machines. If the velocity field is irrotational, the Euler equation (2) 
can be replaced by the Bernoulli equation for the potentials, $„ — (V$i,)^/2 = — 2iJ$„-|-$(,. 
The Z approximation, restricted to the growing mode, requires that each side vanish: one 
side relates $„ and $(, linearly and the other is the "Zel'dovich-Bernoulli" (ZB) equation 
(Nusser & Dekel 1992), 

^,-{D/2){V^,f =0 , (35) 

with the potentials in units of a^ D. The ZB equation can be easily integrated backwards 
with a guaranteed uniform solution at ti. (f^ at to is extractable from observations of 
velocities (§4) or galaxy density (§5), and the initial v and 5 can be derived from the 
initial (fy using linear theory. While the ZB approximation conserves momentum (like 5c) 
one can alternatively satisfy continuity under the Z approximation (like Sq)^ and obtain 
a second-order equation for (fg which can be more accurate than ZB for certain purposes 
(Gramman 1993b, eq. 2.24, 2.25). The recovered initial 5i has deeper valleys and shallower 
hills compared to naive recovery using linear theory, e.g. the GA is less eccentric than 
assumed by Juszkiewicz & Bertschinger (1988). 

Recovering the IPDF. An important issue is whether the initial fluctuations were 
Gaussian or not. A Gaussian fleld is characterized by the joint PDFs of all order being 
generalized Gaussians (cf., Bardeen et al. 1986), and in particular the one-point probability 
of 8 is P{S) oc exp[— (^^/(2(T^)]. Common Inflation predicts Gaussian fluctuations but 
non-Gaussian fluctuations are allowed by certain versions of Inflation (cf., Kofman et al. 
1990) and by models where the perturbations are seeded by cosmic strings, textures, or 
explosions (see Peebles 1993, §16). The present density PDF develops a log-normal shape 
due to non-linear effects (Coles & Jones 1991; Kofman et al. 1994): the tails become 
positively skewed because peaks collapse to large densities while the density in voids cannot 
become negative, and the middle develops negative skewness as density hills contract 
and valleys expand. On the other hand, the PDF of present-day velocity components 
is insensitive to quasi-linear effects (Kofman et al. 1994). 

The observed PDFs today agree with iV-body simulations of Gaussian initial 
conditions (Bouchet et al. 1993), but they have only limited discriminatory power against 
initial non-Gaussianities; the development of a density PDF with a general log-normal 
shape may occur even in certain cases of non-Gaussian initial fluctuations {e.g. Weinberg 
& Cole 1992), and the velocity PDF becomes Gaussian under general conditions due 
to the central limit theorem whenever the velocity is generated by several independent 
density structures. A more effective strategy seems to be to take advantage of the full 
dynamical flelds at to, trace them back in time, and use the linear flelds to discriminate 
between theories. The Eulerian Z approximation can be used to directly recover the initial 
PDF (IPDF) as follows (Nusser & Dekel 1993). The tensor dvi/dxj derived from v{x) 
is transformed to Lagrangian variables q[x). The corresponding eigenvalues /li = dvi/dxi 
and \i = dvi/dqi are related via the key relation 

A, = /i,/(l-/-Vz) • (36) 
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In the Z approximation w oc D so the Lagrangian derivatives A^ are traced back in time 
by simple scahng (xD~^. The initial densities at q{x) can then be computed using linear 
theory, 

<^znOC-(Ai+A2+A3) , (37) 

and the IPDF is computed by bin counting of Sin values across the Eulerian grid, weighted 
by the present densities at the grid points. 

A key feature of the recovered IPDF is that it is sensitive to the assumed value of O 
when the input data is velocities (§8.4), and is 0-independent if the input is density, so that 
it can be used to robustly recover the IPDF from the density field of the 1.2 Jy IRAS survey 
(Nusser et al. 1994). The IPDF turns out to be insensitive to galaxy biasing in the range 
0.5 < 6 < 2, at least for the power-law biasing relation assumed. Errors were evaluated using 
mock IRAS-like catalogs, and the IPDF was found consistent with Gaussian, e.g. the initial 
skewness and kurtosis are limited at the Scr level to —0.65 < S < 0.36 and —0.82 < K < 0.62 
- useful for evaluating specific non-Gaussian models. The first-year COBE measurements 
are consistent with Gaussian but the noise limits their discriminatory power to strongly 
non-Gaussian models (Smoot et al. 1994). 

Gaussianization. A simple method for recovering the initial fluctuations from the 
galaxy distribution under the assumption that they were Gaussian is based on the assertion 
that both gravitational evolution and biased galaxy formation tend to preserve the rank 
order of density in cells. The non-Gaussian distribution of galaxies in cells simply needs to 
be "Gaussianized" in a rank-preserving way (Weinberg 1991). The initial conditions can 
then be evolved forward using an N-body code and compared with the observed galaxy 
distribution in z-space until convergence. A self-consistent solution for a redshift survey in 
the PP region was found if bg ~2, while it was impossible to match the structure on small 
and large scales simultaneously with GI, Gaussian fluctuations and no biasing. This result 
derived independently of O implies a high O once (Sq is determined by another method 
(§8), but note that the limited surveyed region may be an "unfair" sample. 

Least action. The general GI problem with mixed boundary conditions lends itself 
naturally to an application of Hamilton's action principle (Peebles 1989; 1990; 1993). The 
comoving orbit Xi[t) of each mass point nii is parametrized in a way that satisfles the 
boundary conditions 

x,{to)=x,o, lima^i, = 0, (38) 



and the action 



to 




S= I Ldt= I c/t^[(l/2)m,a^i^ -m,$g(s,)] (39) 

i 



is minimized to determine the free parameters. The orbits can be 

X,{t) = X,{to) + Y, fnit)C^,n , (40) 

n 

where fn{t) each satisfy the boundary conditions, and Ci^n ^-re the parameters. The 
problem can be solved to any desired accuracy by increasing n, and a proper choice of 
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/„'s helps the series to converge rapidly to the desired solution. A generalization of the Z 
approximation of the sort 

fn{t) = [D{t) - D{to)r (41) 

is particularly efficient (Giavalisco et al. 1993). In a preliminary application of this scheme 
to a redshift sample (Shaya et al. 1994) the galaxies are assumed to trace mass and be 
self-gravitating. The complication caused by observing redshifts is solved by an iterative 
procedure: a tentative guess is made for the Xi at to based on TF distances to a sub-sample 
of galaxies and a crude flow model, and the least-action solution provides peculiar velocities, 
i.e. redshifts, whose deviations from the observed redshifts are used to correct the Xi for 
the next iteration, until convergence. After the original study of the history of the Local 
Group (Peebles 1989), the method was applied to a sample based on 500 groups within 
30 h~"^Mpc from the Nearby Galaxy Catalog. The results indicate that the unknown tidal 
forces from outside the sampled volume have a signiflcant effect on the recovered flelds 
and should be incorporated in future applications. The solution obtained depends on the 
assumed O, so a comparison with independent TF distances can in principle constrain O 
(§8.2). 

8. THE VALUE OF O 

Assuming that the inferred motions are real and generated by GI, they can be used 
to estimate O. Evidence from virialized systems on smaller scales suggest a low-density 
universe of O ~ 0.1 — 0.2, but these values may be biased. The spatial variations of the 
large-scale velocity fleld now allow measuring the mass density in a volume closer to a 
"fair" sample. One family of methods is based on comparing the dynamical flelds derived 
from velocities to the flelds derived from galaxy redshifts (§9.1, 9.2). These methods can be 
applied in the linear regime but they always rely on the assumed biasing relation between 
galaxies and mass often parametrized by 6, so they provide an estimate of /3 = /(0)/6. 
Another family of methods measures fi from redshift surveys alone, based on 2;-space 
deviations from isotropy (§9.3). Finally, there are methods which rely on non-linear effects 
in the velocity data alone, and provide estimates of O independent of h (§9.4, 9.5). Note 
that the errors quoted by different authors reflect different degrees of sophistication in the 
error analysis, and are in many cases underestimates of the true uncertainty. 

8.1 13 from Galaxies versus the CMB Dipole 

Eq. (28) is best at estimating t>(0), the linear velocity of the LG in the CMB frame 
due to the gravitational acceleration flf(O) exerted by the mass fluctuations around it. A 
comparison with the LG velocity of 627 it 22 kms~"^ as given by the CMB dipole is a direct 
measure of fi. One way to estimate g{Q) is from a whole-sky galaxy survey where only the 
angular positions and the fluxes (or diameters) are observed, exploiting the coincidence of 
nature that both the apparent flux and the gravitational force vary as r~^. If LocM, then 
the vector sum of the fluxes in a volume-limited sample is (xgiQ) due to the mass in that 
volume. This idea can be modifled to deal with a flux-limited sample once the luminosity 
function is known, and applications to the combined UGC/ESO diameter-limited catalog 
of optical galaxies yield fio values in the range 0.3 — 0.5 (Lahav 1987; Lynden-Bell et al. 
1989). These estimates suffer from limited sky coverage, uncertain corrections for Galactic 
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extinction, and different selection procedures defining the north and south samples. The 
IRAS catalog provides a superior sky coverage of 96% of the sky, with negligible Galactic 
extinction and with fluxes observed by one telescope, but with possible under-sampling of 
cluster cores (Kaiser & Lahav 1989). A typical estimate from the angular IRAS catalog is 
l3j = 0.9 ±0.2 (Yahil et al. 1986). 

The redshift surveys provide the third dimension which could help deriving giO) by 
eq. (28), subject to the difficulties associated with discrete, flux-limited sampling (§5). 
The question is whether giO) is indeed predominantly due to the mass within the volume 
sampled, i.e. whether giO) as computed from successive concentric spheres converges 
interior to R^ax- This is an issue of fundamental uncertainty {e.g. Lahav et al. 1990; 
Juszkiewicz et al. 1990; Strauss et al. 1992b). The r—z mapping (21) could either compress 
or rarify the 2;-space volume elements depending on the sign of u in the sense that an outflow 
makes the 2;-space density 8z smaller than the true density 8: 

8z{x) sa 8{x) - 2[v{x) - v{0)] ■ tjr . (42) 

The varying selection function adds to this geometrical effect [in analogy to n{r) in the IM 
bias (§3.2)] and there is contribution from dv I dr as well (Kaiser 1987). It is thus clear that 
the redshifts must be corrected to distances and that any uncertainty in v{x) at large x 
or at ^ = would confuse the derived ^'(0). The latter is the Kaiser "rocket effect": if t?(0) 
originates from a flnite volume r < Vg and the density outside Tq is uniform with w = 0, 
then the measurements in z-space introduce a fake g[0) in the direction of w(0) due to 
the matter outside Tq, and this ^'(0) is logarithmically diverging with r. t?(0) is uncertain 
because it is derived like the rest of v[x) from the density distribution - not from the CMB 
dipole. These difficulties in identifying convergence limit the effectiveness of this method 
in determining the PS on large scales and (3. The hopes for improvement by increasing 
the depth are not high because the signal according to conventional PS models drops with 
distance faster than the shot-noise. 

Attempting to measure (3j from the IRAS data, Strauss et al. (1992) computed the 
probability distribution of g[0) under several models for the statistics of fluctuations, via 
a self-consistent solution for the velocities and an ad hoc flx to the rocket effect, which 
enabled partial corrections for shot-noise, flnite volume, and small-scale non-linear effects. 
They conflrmed that the direction of g[0) converges to a direction only ~ 20° away from 
the CMB dipole, but were unable to determine unambiguously whether \g{0)\ converges 
even within 100 h~"^Mpc. A maximum likelihood flt and careful error analysis constrained 
(3j to the range 0.4 — 0.85 with little sensitivity to the PS assumed. Rowan- Robinson et 
al. (1991, 1993) obtained from the QDOT dipole (3j = O.Sto^^^. Hudson's (1993b) best 
estimate from the optical dipole is fSg =0.72lQ'-^g. 

The volume-limited Abell/ACO catalog of clusters with redshifts within 300 h~"^Mpc 
was used to compute g[0) in a similar way under the assumption that clusters trace mass 
linearly (Scaramella et al. 1991). An apparent convergence was found by ~ 180 h~"^Mpc 
to the value g{0) ~ 4860/3^ km s~^ . A comparison with the LG-CMB motion of 600 km s~^ 
yields (3c ~ 0.123, which corresponds to (3o ~ 0.44 and (3j sa 0.56 if the ratios of biasing 
factors are 4.5 : 1.3 : 1 (§7.1). A similar analysis by Plionis & Valdarnini (1991) yielded 
convergence by ~ 150 h~"^Mpc and (3 values larger by ~30 — 80%. 
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8.2 f3 from Galaxy Density versus Velocities 

The linear correlation found between mass density and galaxy density (§6.2) can be 
used to estimate the ratio (3. The density 5^ determined by POTENT from velocities 
assuming = 1 relates in linear theory to the true 5 hy 5y oc f[Q)5^ while linear biasing 
assumes 5 = b~^5g, so Sy= j38g. Dekel et al. (1993) carried out a careful likelihood analysis 
using the POTENT mass density from the Mark II velocity data and the density of IRAS 
1.9 Jy galaxies, and found [ij = I.SIq'q at 95% confidence. A similar analysis based 
on the Mark III and IRAS 1.2 Jy is in progress. The degeneracy of O and h is broken 
in the quasi-linear regime, where 8{v) is no longer oc f~^. The compatible quasi-linear 
corrections in POTENT and in the IRAS analysis allow a preliminary attempt to separate 
these parameters, which yields for Mark II data O > 0.46 (95% level) if hj > 0.5. A correction 
for IM bias could reduce the 95% confidence limit to O > 0.3 at most. These results are 
valid for linear biasing; possible non-linear biasing may complicate the analysis because it 
is hard to distinguish from non-linear gravitational effects. 

The advantage of comparing densities is that they are local^ independent of reference 
frame, and can be reasonably corrected for non-linear effects. The comparison can 
alternatively be done between the observed velocities and those predicted from a redshift 
survey, subject to limited knowledge of the quadrupole and higher moments of the mass 
distribution outside the surveyed volume and other biases. Kaiser et al. (1991) obtained 
from Mark II velocities versus QDOT predictions f3j = 0.9lo'i5- An analysis by Roth 
(1994) using IRAS 1.9 Jy galaxies yielded (3j = 0.6 ± 0.3 {2a). Nusser & Davis (1994) 
implemented a novel method based on the Z approximation in spherical harmonics to 
derive the velocity dipole of distant shells from the IRAS 1.2 Jy survey and found in 
comparison to the dipole of observed velocities (3j = 0.6 it 0.2. 

Similar comparisons with the optical galaxy fields indicate a similar correlation 
between light and mass. A comparison at the velocity level gives (Hudson 1994) (Sq = 
0.5 lb 0.1, and a preliminary comparison at the density level with 12 h~"^Mpc smoothing 
indicates (Hudson et al. 1994) (3o ~ 0.6 it 0.2, in general agreement with the ratio of 
bo /bj sa 1.3—1.4 obtained by direct comparison. Shaya et al. (1994) applied the least-action 
reconstruction method (§7.2) to a redshift survey of several hundred spirals within our local 
30 h~"^Mpc neighborhood and crudely obtained by comparison to TF distances (3o ~0.4. 

8.3 f3 from Distortions in Redshift Space 

Redshift samples, which contain hidden information about velocities, can be used on 
their own to measure (3. The clustering, assumed isotropic in real space, x, is anisotropic in 
z-space^ z, where z = r-\-x-'V displaces galaxies along the preferred direction X. While virial 
velocities on small scales stretch clusters into "fingers of god" along the LOS, systematic 
infall motions enhance large-scale structures by artificially squashing them along the LOS. 
The linear approximation —^■v = (35g indicates that the effect is (3 dependent because — V-« 
is related to the anisotropy in z-space while 5g is isotropic, so the statistical deviations 
from isotropy can tell (3 {e.g. Sargent & Turner 1977). 

Kaiser (1987) showed in linear theory that the anisotropic Fourier PS in z-space is 
related to the r-space PS of mass density, P{k), via 

P^{k,^,) = P{k){l + f3^,y , (43) 
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where ii = k-x. This relation is vahd only for a fixed /i, i.e. in a distant volume of small 
solid angle (Zaroubi & Hoffman 1994), but there are ways to apply it more generally. The 
redshift PS can be decomposed into Legendre polynomials, Vi{ii)^ with even multipole 
moments Pf[k)^ 

P^{k,^) = Y^Pl{k)Vi{p:) , Pl{k) = ^- / d^P^{k,^)Vi{ii) . (44) 

Based on eq. (43) the first two non- vanishing moments are 

P^{k) = (1 + ^/3 + \f)P{k) , P^{k) = {^(3 + ^(3')P{k) , (45) 

so the observable ratio of quadrupole to monopole is a function of (3 independent of P{k). 
A preliminary application to the 1.2 Jy IRAS survey yields (3j ~ 0.3 — 0.4 at wavelength 
30 — 40 h~"^Mpc, suspected of being an underestimate because of non-linear effects out 
to ~ 50 h~"^Mpc (Cole et al. 1993). Peacock & Dodds (1994) developed a method for 
reconstructing the linear PS and they obtain /3j = 1.0 it 0.2. 

The distortions should be apparent in the 2;-space two-point correlation function, 
(^^(rp,7r), the excess of pairs with separation tt along the LOS and r^ transversely (Davis 
& Peebles 1983). The contours of equal (^, assumed round in r-space, appear in 2;-space 
elongated along the LOS at small separations and squashed on large scales depending on 
fi. Hamilton (1992; 1993) used the multiple moments of (^^, in analogy to eqs. (26-27), 
and his various estimates from the 1.9 Jy IRAS survey span the range fij = 0.25—1. Fisher 
et al. (1994a) computed C^(rp,7r) from the 1.2 Jy IRAS survey, and derived the first two 
pair-velocity moments. Their attempt to use the velocity dispersion via the Cosmic Virial 
Theorem led to the conclusion that this is a bad method for estimating O, but the mean, 
(^12) = 109147 ^^ 10 h~"^Mpc, yielded /3j = 0.45lo'i8- The drawbacks of using ^ versus PS 
are that (a) the uncertainty in the mean density affects all scales in ^ whereas it is limited 
to the A; = mode of the PS, (b) the errors on different scales in ^ are correlated whereas 
they are independent in a linear PS for a Gaussian field, and (c) ^ mixes different physical 
scales, complicating the transition between the linear and non-linear regimes. Non-linear 
effects tend to make all the above results underestimates. 

A promising method that is tailored to deal with a realistic redshift survey of a 
selection function </)(r) and does not rely on the subtleties of eq. (43) is based on a 
weighted spherical harmonic decomposition of (^^(z) (Fisher et al. 1994b), 

at^ = jd'z^{r)f{z)[l+5,{z)]YUs). {\aU') = lJ^ dkk'P{k)\4,l{k) + (34,f{k)\' . 

(46) 
The arbitrary weighting function f[z) is vanishing at infinity to eliminate surface terms. 
The mean-square of the harmonics is derived in linear theory assuming that the survey is 
a "fair" sample, and rp^ and rp'^ are explicit integrals over r of certain expressions involving 
(/)(r), /(r), Bessel functions and their derivatives. The first term represents real structure 
and the second is the correction embodying the 2;-space distortions. The harmonic PS in 
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z-space^ averaged over m^ is thus determined by P{k) and /3, where the 2;-space distortions 
appear as a /3-dependent excess at smaU /. The harmonic PS derived from the 1.2 Jy IRAS 
survey yield (3j = 1.0 it 0.3 for assumed erg = 0.7 (motivated by the IRAS (^, Fisher et al. 
1994a), with an additional systematic uncertainty of ±0.2 arising from the unknown shape 
of the PS. 

The methods for measuring fi from redshift distortions are promising because they are 
relatively free of systematic errors and because very large redshift surveys are achievable 
in the near future. With a sufficiently large redshift survey, one can even hope to be able 
to use the non-linear effects to determine O and h separately. 



8.4 O from PDFs using velocities 

Assuming that the initial fluctuations are a random Gaussian fleld, the one-point PDF 
of smoothed density develops a characteristic skewness due to non-linear effects early in 
the quasi-linear regime (§7.2). The skewness of 8 is given in second-order perturbation 
theory by 

{8^)l{8y = {MlT-?,-n) , (47) 

with n the effective power index near the smoothing scale (Bouchet et al. 1992). Since this 
ratio for 8 is practically independent of O, and since V-«~ —f8, the corresponding ratio 
for V-t> strongly depends on O, and in second-order (Bernardeau et al. 1994) 

S, = {{V-vf)/{{V-vf)' = -/(0)-i(26/7 - 3 - n) . (48) 

Using N-body simulations and 12 h~"^Mpc smoothing one indeed flnds Ss = —1.8 it 0.7 for 
0=1 and 5*3 = —4.1 it 1.3 for O = 0.3, where the quoted error is the cosmic scatter for a 
sphere of radius 40 h~"^Mpc in a CDM universe [Hq = 75, 6=1). A preliminary estimate 
of Ss in the current POTENT velocity fleld within 40 h~"^Mpc is —1.1 it 0.8, the error 
representing distance errors. With the two errors added in quadrature, = 0.3 is rejected 
at the ~2(T level (somewhat sensitive to the assumed PS). 

Since the PDF contains only part of the information stored in the data and is in 
some cases not that sensitive to the IPDF (§7.2), a more powerful bound can be obtained 
by using the detailed v[x) to recover the IPDF, and use the latter to constrain O. This 
is done by comparing the 0-dependent IPDF recovered from observed velocities to an 
assumed IPDF(Nusser & Dekel 1993), most naturally a Gaussian as recovered from IRAS 
density (§7.2). The velocity out of POTENT Mark II within a conservatively selected 
volume was fed into the IPDF recovery procedure with O either 1 or 0.3, and the errors 
due to distance errors and cosmic scatter were estimated. The IPDF recovered with = 1 
is found marginally consistent with Gaussian while the one recovered with O = 0.3 shows 
signiflcant deviations. The largest deviation bin by bin in the IPDF is ~2(T for = 1 and 
>4(T for O = 0.3, and a similar rejection is obtained with a X^"type statistic. The skewness 
and kurtosis are poorly determined because of noisy tails but the replacements (a;|a;|) and 
{\x\) allow a rejection of O = 0.3 at the (5 — 6)0" levels. 

37 



8.5 O from Velocities in Voids 

A diverging flow in an extended low-density region can provide a robust dynamical 
lower bound on O, based on the fact that large outflows are not expected in a low-0 
universe (Dekel & Rees 1994). The velocities are assumed to be induced by GI, but no 
assumptions need to be made regarding galaxy biasing or the exact statistical nature of the 
fluctuations. The derivatives of a diverging velocity fleld infer a non-linear approximation 
to the mass density, 8(.{^,dv I dx) (eq. 7), which is an overestimate, 8(.>8, when the true 
value of O is assumed. Analogously to 8q = — /(0)~"^V-w, the S^. inferred from a given 
diverging velocity fleld becomes more negative when a smaller O is assumed, and it may 
become smaller than —1. Since 8 > —1 because mass is never negative, O is bounded from 
below. 
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Figure 8: Maps of 8c inferred from the observed velocities near the Sculptor void 
in the Supergalactic plane, for two values of O. The LG is marked by '-|-' and the 
void is conflned by the Pavo part of the GA (left) and the Aquarius extension of 
PP (right). Contour spacing is 0.5, with 8,. = Q heavy, 8,. > Q solid, and 8,. < Q 
dotted. The heavy-dashed contours mark the illegitimate downward deviation of 
8c below —1 in units of cr^, starting from zero [i.e. 8c = —1), and decreasing with 
spacing — O.Scr. The value O = 0.2 is ruled out at the 2.9(T level (Dekel & Rees 
1994). 

The inferred 8c{x) smoothed at 12 h~"^Mpc and the associated error fleld a§ are derived 
by POTENT from the observed radial velocities and, focusing on the deepest density 
wells, the assumed O is lowered until 8c becomes signiflcantly smaller than —1. The most 
promising "test case" provided by the Mark III data seems to be a broad diverging region 
centered near the supergalactic plane at the vicinity of (X, Y) = (—25, —40) in h~"^Mpc - 
the "Sculptor void" of galaxies (Kauffman et al. 1991) next to the "Southern Wall" (Figure 
8). Values of Osa 1 are perfectly consistent with the data, but 8c becomes smaller than —1 
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already for = 0.6. The values = 0.3 and 0.2 are ruled out at the 2.4-, and 2.9cr levels 
in terms of the random error as- This is just a preliminary result. The systematic errors 
have been partially corrected for in POTENT, but a more specific investigation of the SG 
biases affecting the smoothed velocity field in density wells is required. For the method to 
be effective one needs to find a void that is (a) bigger than the correlation length for its 
vicinity to represent the universal O, (b) deep enough for the lower bound to be tight, (c) 
nearby enough for the distance errors to be small, and (d) properly sampled to trace the 
velocity field in its vicinity. 



TABLE 1: ft and i* 



CMB dipole 


VS galaxies angular 


Yahil et al. 86 


A- 


:0.9±0.2i 




VS galaxies redshift 


Strauss et al. 92 


A- 


: 0.4 -0.85 






Rowan-Rob. et al. 91 


A- 


■ 8+°-2 

"-'■°-0.15 




VS galaxies angular 


Lynden-Bell et al. 89 


/3o- 


^0.3-0.5 




vs galaxies redshift 


Hudson 93b 


[io = 


- 7+°-'' 

- '-'■'-0.2 




clusters 


Scaramella et al. 91 


[ia ~0.13 






Plionis et al. 91 


fia ~ 0.17 -0.22 


V VS 5g 


Potent-IRASl.9 density 


Dekel et al. 93 


A- 


: 1.31°:^^ (95%) 




Potent-IRASl.2 v-dipole 


Nusser h Davis 94 


A- 


: 0.6 ±0.2 




TF-QDOT 


Kaiser et al. 91 


A- 


■ 9+°-^ 

'-'■^-0.15 




TF-QDOT clusters 


Frenk et al. 94 


A- 


: 1.0 ±0.3 




TF inverse - IRAS 1.9 


Roth 94 


A- 


: 0.6 ±0.35 (2cr) 




Potent- Optical density 


Hudson et al. 94 


/3o- 


^0.8 ±0.2 




TF-Optical 


Hudson 93c 


/3o- 


^0.5 ±0.1 




TF-Optical local 


Shaya et al. 94 


fio ~ 0.4 


5;-distortions 


Pk IRAS 1.2 


Cole et al. 93 


A> 


,0.3-0.4 




Pk IRAS 1.2 


Peacock L Dodds 94 


A- 


: 1.0 ±0.2 




C IRAS1.9 


Hamilton 93 


A- 


'0.25- 1 




^ IRAS1.2 


Fisher et al. 94a 


A- 


: 0.451°:^ 




Yim IRAS 1.2 


Fisher et al. 94b 


A- 


: 1.0 ±0.3 


Velocities 


Gaussian IPDF Potent 


Nusser & Dekel 93 


VL> 


0.3 (4ct - 6ct) 




Skew(V-«) Potent 


Bernardeau et al. 94 


Q> 


0.3 (2ct) 




Voids Potent 


Dekel & Rees 94 


Q> 


0.3 (2.4ct) 



l3 = Q"-^/b, be -.bo-.bj ^ 4.5 : 1.3 : 1.0, see text. 
All errors are la unless stated otherwise. 
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9. DISCUSSION: ARE THE HYPOTHESES JUSTIFIED? 

Are the motions real? The large-scale bulk flow rests upon the interpretation of 
the CMB dipole as motion of the LG in the CMB frame, based on the assumption of 
cosmological isotropy (HI), and supported by the gravitational acceleration derived at the 
LG from the galaxy and mass distribution around it (§6.1). The detection of 5T/T ^ 10~^ 
is a clear indication, via GI, for ~ 300 kms~"^ motions over ~ 100 h~"^Mpc (§6.1). While 
this is reassuring, the LP hint for possible very-large coherence suggests caution (§7.1). 
The evidence in support of H6, that the TF-inferred motions about the LG are real, are 
(a) the correlation 5g oc — V-t>, which is robustly predicted for true velocities based on 
continuity and would be hard to mimic by environmental effects (§6.2), (b) the failure to 
detect any signiflcant correlation between velocities and the environment or other galaxy 
properties (§6.3), and (c) the similarity between the velocity flelds traced by spirals (TF) 
and by ellipticals (D„— cr) (§6.3). 

Is linear biasing a good approximation (H5)? The galaxy- velocity correlation is most 
sensitive to it, and the observed correlation on scales ^10 h~"^Mpc is consistent with 
linear biasing properly modifled in the tails (§6.2). However, it is diflicult to distinguish 
non-linear biasing from non-linear gravitational effects, and the range of different estimates 
of (3 (§8) may indicate that the biasing parameter varies as a function of scale. The ratio 
of ~ 10 h~"^Mpc-smoothed densities for optical and IRAS galaxies is bg /bj sa 1.3 — 1.5. 

Is gravity the dominant source of LSS (H4)? The observed velocity-density correlation 
(§6.2) is fully consistent with GI, but it is sensitive to continuity more than to the speciflc 
time dependence implied by gravity. Any non-GI process followed by a gravitating phase 
would end up consistent with this observation, and certain non-GI models may show a 
similar spatial behavior even if gravity never plays any role. The E-S correlation (§6.3) is 
also consistent with gravity as galaxies of all types trace the same velocity fleld (H4b), 
but any model where all galaxies are set into motion by the same mechanism could 
pass this test. The strongest evidence for gravitational origin comes from the statistical 
agreement between the fluctuations of today and those implied by the CMB at the time 
of recombination. A marginal warning signal for GI is provided by the ~ 700 kms~"^ 
bulk velocity indicated by LP for rich clusters across ~ 200 h~"^Mpc. Such a velocity at 
face value would be in conflict with the gravitational acceleration implied by the cluster 
distribution and with the 5T/T^ 10~^ at ~2°, but the errors are big. 

The property of irrotationality (H4a) used in the reconstruction from either velocities 
or densities is impossible to deduce solely from observations of velocities along the lines of 
sight from one origin. Irrotationality is assumed based on the theory of GI, or it can be 
tested against the assumption of isotropy by measuring the isotropy of the velocity fleld 
derived by potential analysis for a fair sample. 

The observed CMB fluctuations provide evidence for initial fluctuations (H2), 
consistent with a scale-invariant n ~ 1 spectrum. The observed motions are also consistent 
with n~ 1 (§7.1) (with the uncertain LP result as a possible exception). The indications 
for somewhat higher large-scale power in the clustering of galaxies (Maddox et al. 1990) 
may reflect non-trivial biasing. The question of whether the fluctuations were Gaussian 
(H2) cannot be answered by observed velocities alone. The PDF of V • w is consistent 
with Gaussian initial fluctuations skewed by non-linear gravity, but this is not a very 
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discriminatory test. Nevertheless, the galaxy spatial distribution does indicate Gaussian 
initial fluctuations fairly convincingly (§7.2) 

Can we tell the nature of the dark matter (H3)? In view of the tight nucleosynthesis 
constraints on baryonic density, the high O indicated by the motions requires non-baryonic 
DM. The mass-density PS on scales 10 — 100 h~"^Mpc is calculable in principle but the 
current uncertainties do not allow a clear distinction between the possibilities of baryonic, 
cold, hot or mixed DM. The mixed model seems to score best in view of the overall LSS 
data, as expected from a model with more free parameters, but CDM in fact does somewhat 
better in fltting the large-scale motions. I do not think that any of the front-runner models 
is signiflcantly ruled out at this point, contrary to occasional premature statements in the 
literature about the "death" of certain models. I predict that were the DM constituent (s) 
to be securely detected in the lab, the corresponding scenario of LSS would flnd a way to 
overcome the ~ 2cr obstacles it is facing now. 

What can we conclude about the background cosmology (HI)? All the observations so 
far are consistent with large-scale homogeneity and isotropy (with the exception of the 2cr 
LP discrepancy). The motions say nothing about H or A (Lahav et al. 1991), but they 
provide a unique opportunity to constrain O in several different ways. Some methods put a 
strong (>3(t) lower bound of 0> 0.2—0.3. This is consistent with the theoretically-favored 
= 1 but "ugly" values near O sa 0.5 are not ruled out either. The apparent range of fi 
values obtained on different scales is partly due to errors, and any remaining difference may 
be explained by a scale-dependent non-linear biasing relation between the different galaxy 
types and mass. The data are thus at least consistent with the predictions of Inflation: 
flat geometry and Gaussian, scale-invariant initial fluctuations. Recall however that = 1 
predicts to = 6.3/z~"^Gyr, which will be in conflict with the age constraints from globular 
clusters, to = 15 it 3Gyr, if the Hubble constant h is not close to 0.5. 

The rapid progress in this fleld guarantees that many of the results and uncertainties 
discussed above will soon become obsolete, but I hope that the discussion of concepts will 
be of lasting value, and that the methods discussed can be useful as are and as a basis for 
improvements. 
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